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ABSTRACT 

Two observed features of dwarf spheroidal (dSph) galaxies — extended light profiles and more or less flat 
velocity dispersion profiles — have led to competing interpretations over the nature of dSphs: Either the outer 
characteristics of dSphs are reflecting ongoing tidal disruption, or they are showing the existence of bound stel- 
lar "halos" embedded within extended dark matter halos. To accommodate the observed dSph features within 
the latter interpretation, dSph models now commonly adopt multi-component structures with strongly varying 
M/L with radius. In contrast, we explore whether the observed features of dSphs can still be accounted for with 
simpler, mass-follows-light (MFL) models including tidal disruption. As a test case, we focus on modeling the 
Carina dSph, which resembles other classical dSphs in shape and velocity distribution but presently has the 
most extensive data at large radius. We find that previous A/-body, MFL simulations of (tidally disrupting) 
dSphs did not sufficiently explore the parameter space of satellite mass, density and orbital shape to find ade- 
quate matches to Galactic dSph systems, whereas with a systematic survey of parameter space we are able to 
find tidally disrupting, MFL satellite models that rather faithfully reproduce Carina's velocity profile, velocity 
dispersion profile and projected density distribution over its entire sampled radius. The successful MFL model 
satellites have very eccentric orbits, currently favored by Cold Dark Matter (CDM) models, and central veloc- 
ity dispersions that still yield an accurate representation of the bound mass and observed central M/L ~ 40 
of Carina, despite inflation of the velocity dispersion outside the dSph core by unbound debris. Our survey 
of parameter space also allows us to address a number of commonly held misperceptions of tidal disruption 
and its observable effects on dSph structure and dynamics. The simulations suggest that (1) even modest tidal 
disruption can have a profound effect on the observed dynamics of dSph stars at large radii, and (2) it may be 
premature to abandon single-component, MFL models, which can successfully reproduce all salient structural 
and dynamical properties of the Carina (and other) dSph satellites without invoking special timing or viewing 
perspectives. Satellites that are well-described by tidally disrupting MFL models could still be fully compatible 
with A-CDM if for example they represent a later stage in the evolution of luminous subhalos. 

Subject headings: Carina Dwarf - galaxies: stellar content - Local Group - kinematics and dynamics - 
cosmology : dark matter 



1. INTRODUCTION 

1.1. Dark Matter and Tidal Effects in dSphs: Previous Work 

Ever since the work of Hod ge and Michie dHodgel 1961 alibi 
Il962t iHodge & Michiel lT969). which suggested that Galac- 
tic tides might play an important role in shaping the struc- 
ture of Galactic dwarf spheroidal galaxies (dSphs), the ex- 
tent of this tidal influence has been debated. When early 
precision radial velocity measurements wer e obtained in the 
central regions of a number of these dSphs dAaronsonll 19831 : 
iMateo et~afl 119931: lHargreaves et al.1 1 19941) the relative im- 
portance of tidal forces became less clear because the mea- 
sured high velocity dispersions implied large dark matter con- 
tents within dSphs. Using traditional core-fi tting methods 
(e.g.. Ulling worth O 976t iRichstone & Tremainel l 19861) . central 
mass-to-light ratio (M/L) values ra nging from ~ 5 for Leo I 
( IMateo et alJI 1 9981: ISohn et al.ll2007h to - 100 for Draco and 
Ursa Minor (Armandroff et al. 1995), and even a s high as 500 
for Ursa Major ( Klevna et al. 2005) and Bootes dMunoz et all 
2006b; iMartin et alj|2007l) have been found among Galactic 
dSph galaxies. 

Such large M/L's have been difficult to reconcile with 
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the notion that these dwarfs may be tidally truncated 
(ICuddeford & Millerffl990l lBurkertlll997h . Moreover, were 
they to have such high M/L, dSphs would be structurally 
very different from systems of relatively similar mass and 
luminosity scale, such as globular clusters and dwarf ellip- 
tical galaxies. The notion that dSph velocity dispersions 
may not directly reflect their mass but may be inflated by 
various mechanisms has prompted frequent assessment of 
various alternative dSph structural models — both New- 
tonian and n on-Newtonian (e.g., MOND; Milgrom 11995 . 
i Lokasi l2001l but see also, e.g.. | Gerhard & Spergel 11992 , 
Sanchez-Salcedo & Hernande 3 12001 " — as well as the as- 
sumptions built into the core-fitting method for deriving 
M/L's. 4 

For example, among the suggestions within a Newtonian 
framework, the possible influence of binary stars to inflate 
velocity dispersions was investigated, but found to be minor 
dOlszewski et all 1 19961: lHargreaves etalj|1996t iKlevna et all 
|2002|) . The influence of velocity jitter in the atmospheres 
of giant stars is also likely i nsignificant since suc h variations 
rarely exceed .5-1 km s" 1 dBizyaev et al.ll2006l) or perhaps 
slightly more dCarnev et al.l 12003k Other explanations for 
the seemingly high observed dSph velocity dispersions ques- 

4 Core-fitting assumes that these dSphs are in virial equilibrium, their ve- 
locity distributions are isotropic, their mass distribution follows that of the 
light, and their surface brightness, or equivalently, their surface density, are 
well represented by a King 1 1962) model. 
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tion the validity of the virial equilibrium assumption, in some 
cases exploring violations of that equilibrium in the extreme 
— i.e. proposed models whereby the dSphs are in com- 
plete or near complete tidal disruption and their appearance of 
structural cohesion is an artefact of other dynamical processes 
that obviate the need for dark matter. Synchronicity between 
internal stellar and external bulk dSph orbits, for instance, 
has been proposed as one means to explain the high veloc- 
ity dispersions and apparent spatial coheren ce of dSph sys- 
tems i n the resonant orb it coupling models of Kuhn & Miller 
( 1989) andlKuhnldl993h. Sub sequent work by iPrvorf jl 996ft 
and Sellw ood & Pryorl (1 1998) showed that the tidal oscilla- 
tion mechanism proposed by lKuhn & Milled (11989ft is not ex- 
cited by motion through a realistic, logarithmic Galactic po- 
tential, andjhat the tidally disrupted remnants proposed by 
iKuhnl (11993ft would be some ten times large r than the ob- 
ser ved dSphs. However, Fleck & Kuhnl d2003l) have updated 
the iKuhn & Miller! dl989h model and suggest that tidally in- 
duced oscillations can inflate the velocity dispersion by up to 
an order of magnitude, although only for a short fraction of a 
dSph orbital period. 

Meanwhile, models of dSphs that are mostly bound 
but that incorporate the presence of some unbound stars 
have been partly motivated by the discovery that most 
of the Galactic dSphs have an extended component be- 
yond their nominal King limiting radius r;,,„ (obtained 
by fitting a King profile to the inner light distributions; 
see Appendix A), a feature suggeste d to be eviden ce 
of tidal disruption in these systems (Eskridg g Il988alfbt 
Irwin & Ha tzidimitri ouiri99l IKuhn. Smith & Hawleyl [1991 
Smith. Kuhn, & Hawlevlll997t iMaiewski et al.ll2000t) Since 
these early photometric studies, extended density profile 
"break" populati ons 5 have been found and/o r verified in 
the U r sa Minor dMartfnez-Delgado et al.ll2001l lPalmaetal.1 
20031: iMufioz et al.l 120051). Sculptor dWalcher et al.1 120031: 
Westfall et al.1 \20od . Draco dWilkinson et al.l 120041). Leo I 
dSohn et alJl2007l) and Car ina dSphs dMajewski et al.1 120001 
2005; Munoz et al. 2006a). As has been pointed out in 
many of these references, the observed "bimodal" dSph 
structural profiles in these systems look very similar to 
those produced by jV-b ody simulations of tidally disrupt- 
ing dwarf galaxies (e.g.. Ijohnston et al.l 1 1 999bl iMaver et ail 
120021: IChoi. Guhathakurta. & Johnstonll2002l) . and in at least 
one other real dSph case s haring this radial profi le shape — 
the Sagittarius (Sgr) dSph (IMaiewski et al.ll2003l) — the sec- 
ondary, extended population is unquestionably constituted by 
tidally stripped stars. By analogy it might seem logical to 
surmise that the other dSph profiles have been shaped by the 
same disruption processes. However, two frequently cited N- 
body studies of dSph tidal disruption are typically invoked to 
argue against this conclusion. 

In an attempt to investigate the degree to which Galactic 
tides can affect the structural and kinematic al properties of 
dSphs - and therefore their inferred M/L - Piatek & Pryoj 
(fl99l hereafter PP95) and lOh. Lin. & Aarsethl (fl99l here- 
after OLA95) each carried out numerical simulations of low 
density satellites containing no dark matter orbiting in a Milky 
Way (MW)-like potential. IOLA95I modeled low mass satel- 
lites (M sate iiite < 2 x 10 6 M ) orbiting in rather dynamically 

5 Here, "break" population refers to an empirical description of the density 
profile where the outer parts of the density distribution seems to "break" to a 
new, less steep form, departing from a distribution (usually resembling a King 
profile — see Appendix A) that fits well the inner parts of Galactic dSphs. 



gentle orbits (eccentricities no greater than 0.5, and perigalac- 
tica no closer than 50 kpc) having several perigalacticon pas- 
sages. OLA95 found that for such systems the satellite den- 
sity profile becomes flattened in the orbital plane of the galaxy 
and follows aproximately a power-law distribution with long 
strand s of tida lly stripped debris formed in some cases. How- 
ever, IOLA95I also found that the velocity dispersion of the 
unbound, but not yet dispersed, stars was similar to the pre- 
vious central equilibrium value, implying that the observed 
high velocity dispersions of dSph s were not being signifi- 
cantly inflated by tidal stripping; IOLA95I concluded, there- 
fore, that substantial mass content within the dSphs is un- 
avoidable, even when/if tides are important. 

PP95 also studied low mass satellites with masses of 2.7 x 
10 5 M Q and 8.1 x 10 s M Q but in rather eccentric orbits 
(apogalactica of 210 kpc and perigalactica of 30 kpc). Again, 
rat her stro ngly formed tidal tails were formed in these models, 
but lPP95l found that a significant, artificial (i.e., tidal) increase 
of the M/L was rarely observed, and, even when found, the 
phenomenon lasts only briefly (several times 10 8 yr); to in- 
voke this scenario as an explanation for the high observed 
central velocity dispersion in dSphs would require that we 
live in a very special time to be observing all dSphs in the 
same evolutionary phase (near perigalacticon). Thus, IPP95I 
concluded that it is very unlikely that the central regions of 
dSphs are influenced by tides. They also pointed out that the 
effects of Galactic tides on dSphs were, for the most part, un- 
detectable, except perhaps for two key features: (1) a velocity 
gradient (apparent rotation) that should be observed along the 
major axis of tidally disrupting dwarf galaxies and, (2) a sur- 
face brightness distribution that should become more elliptical 
with r adius. 

The IOLA95I and |PP95 studies have been highly influen- 
tial in reinforcing the notion that tidal disruption is not oc- 
curing in most dSph systems, since observed dSph galax- 
ies with large central velocity dispersions do not resemble 
these particular models of disrupted systems, which showed 
tidally induced rotation, distorted morphologies and low cen- 
tral velocity dispersions. Thus has evolved a frequently 
presented, "bimodal" pair of viewpoints that dSph galax- 
ies either have little dark matter and therefore are highly 
prone to disrupting, or that they have high dark matter con- 
tents that are incommensura t e with tidal stripping 6 (e.g. , 
[ Burked 119971; iKroupal 119971; iGom ez-Flec hoso etall 119991: 
Kleynaetal. 2001; Gomez-Flechoso & Martmez-Delgado 
l2Q03t iGilmore et al.1 120061: iKoch et al.l l2007al). Ho wever, 
though not technically flawed, the IOLA95I and IPP95I studies 
are highly specialized — the former pertaining spec ifically to 
low mass dSphs on rather circular orbits whereas [PP95I se- 
lected as strawman examples models with M/L so unrealis- 
tically low that the satellites do not even survive one peri- 
galacticon passage — so that it is not clear how general are 
the conclusions of these studies and whether they should con- 
tinue to be applied to the interpretation of the much better data 
we now have on dSph systems. For example, in reality it is 
likely that dSphs have rather eccentric orbits, as is now known 

6 Even when dark matter and tides are included i n models, t he con- 
clus ions lead to a bimodal position. For example, Read et al. 1 2006b) 
and Mashchenko et al. (2006) combined tidal interaction with dark matter- 
dominated, two-component satellite models when studying Galactic dSphs, 
but concluded that tides are unimportant due to the immense amounts of dark 
matter in the satellite extended dark matter halos they need to explain the 
velocity dispersion profiles of dSphs, effectively reinforcing the "bimodal" 
picture of tides versus dark matter. 
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to be the case f or at least severa l of the most recent orbits of 
the Sgr system (|Law et al.ll2005l) . as is expected for the Leo I 
system (with its enormous rad ial velocity), as has been mea- 
sured for several other dSphs dPiatek et al J20031 12005L 120061) 
and as is predicted f or cold dark matter (CDM) subhalos (e.g., 
iGhigna et ai1ll998l) . Moreover, the PP95-adopted M/L = 1 
Mq/Lq is several times lower than that expected even for an 
old stellar population with no dark matter, so that the IPP95I 
modeled satellites are much fluffier than would be expected 
for any dSph system, even ignoring dark matter. While im- 
portant for demonstrating how dwarf satellites would evolve 
in specific circumstances, these early models should be used 
cautiously as guides to understanding real dSph systems on 
eccentric orbits. 

Subsequent to these 1995 studies, iRroupal ( 1 1997b and 
iKlessen & Kroupal dl998) were successful in explaining many 
of the salient features observed in present-day dSphs by mod- 
eling them as long-lived, dark matter free, but mostly unbound 
remnant dSph cores of mass only a few percent of the ini- 
tial satellite mass. In these simulations, the satellites have 
physical dimensions comparable to Galactic dSphs, they fol- 
low the observed linear c orrelation between sur face bright- 
ness and total luminosity (Bellazzi ni et al.lfl996l) . and more 
importantly, they have inferred (but artificial) M/L as high 
as 100. However, these models only work if the satellites 
are in such an eccentric orbit that we see them and their 
debris more or less projected along the same line of sight. 
Such alignments are directly testable because the projection 
of tidal debris is expected to result in a significant spread in 
color m agnitude diagram (CM P) features like the horizontal 
branch (IKlessen & Zhaoir2002l) . This depth effect should be 
detectable with wide-field photom etry, a test carried out by 
IKlessen. Grebel. & Harbeckl (f2003) with public data from the 
Sloan Digital Sky Survey (SDSS) in the Draco dSph field. 
They showed that Draco's narrow horizontal branch cannot 
be reproduced by this alignment model and concluded that 
Draco cannot be the unbound remnant of a disrupted galaxy, 
and must, therefore, be dark matter dominated. This addi- 
tional "failure of the tidal model" — though again pertain- 
ing to a prediction from a rather highly specialized satellite 
model — has further promoted a perception that tidal effects 
are unimportant in dSphs or at least h a ve not affected them 



significantly (e.g. , iKlevna et al 



iRead et alj 2006b: Gilmo re et al 



l20Qlt iWalker et al.1 l2006a: 



12003) . 



More recently, new, large area radial velocity (RV) datasets 
have become available for the Galactic dSphs and allow for 
the first time the study of the line of sight velocity dispersion 
as a function of radius out to (and in some cases well past) 
r/, m . Such extensive velocity data provide new and remark- 
ably consistent constraints on the structure of dSphs. Flat or 
only slightly rising/declini ng velocity dispersion profiles have 
been reported f or Sculptor dTolstov et al.|l2004UWestfall et all 
120061) . Drac o dWilkinson et alj 12004 iMunoz et al.1 120051) . 
Ursa Minor dWilkinson et al.T l2004: Mu noz et al 112005b For- 
na x (IWalker et al . 2006a), Sextans (IWalker et al.ll2006bT). Leo 



I (Sohnetal. 2007 



2007b, Siegeletal 



Koc h et al.1 l2007al). LeoTT dKoch et all 
120071) . Carina TMunozetal.ll2006al) and 
Sagittarius (S. R . Majewski et al., in preparation ). Previ- 
ous studies (e.g., lKroupa|[l997l: IKlevna et al.ll 19991) have de- 
mostrated that an important difference between a dark matter- 
dominated, non-disrupting, mass-follows-light (MFL) dSph 
and a tidally disrupted one is the radial behavior of the veloc- 
ity dispersion. The former model would yield a decreasing 
dispersion with radius, whereas a flat or rising profile is ex- 



pected for a tidally disrupting dwarf. Thus, tidal disruption 
provides one explanation for the observed dSph dispersion 
profiles. 

But a third structural possibility, a luminous dSph sur- 
rounded by a more extended, mostly dark matter halo is 
also consistent with the new observations. CDM cosm ology- 
driven models dStoehr et al.ll2002t lHavashi et aT]|2003l) have 
shown that the observed dSph velocity dispersion profiles 
are in good agreement with simulations wherein dSphs in- 
habit the most massive subhalos (10 9 - 10 10 M Q ) in a A- 
CDM universe, a result that also helps mitigate the so- 
calle d missing satellite problem inherent to these cosmolo- 
gies dKauffmann et al.lfT99l iKlypin et ai1ll999t iMobre et all 
1999). This theoretical result, and the apparent lack of ob- 
servational evidence supporting tidal str ipping a s an ongoing 
process (in the context of the PP95 and OLA95 models), has 
prompted some more recent studies (e. g .. IKl evna et al. ]|2002 ; 
Lokasll200l iM ashchenkoeTall l2006t IWalker et al.1 12006a: 



Read et al . 2006b) to model the structure of dSphs with two 



component models (compact stellar and extended dark mass 
components) to account for the observed dSph velocity dis- 
persion profiles; these two component models use extended 
DM halos to keep the velocity dispersion profiles flat, and then 
conveniently create total masses consistent with the enormous 
total M/L [exceeding 10 3 or 10 4 (M/L)©] required for dSphs 
to populate the high end of the subhalo mass spectrum in A- 
CDM predictions. 7 These models explicitly assume that the 
dwarf satellite is in virial equilibrium thoughout its physical 
extent, and therefore, neglect any tidal influence in the shape 
of the velocity distribution. 

1.2. Goals of this Paper 

Against this backdrop of previous work on dSph structure 
and tidal influences, the present paper has several goals. First, 
we want to readdress the question to what extent Galactic 
tides may be affecting dSph galaxies, and, more specifically, 
reassess the issue of whether other presently observed Galac- 
tic dSph satellites may be tidally disrupting along the Sgr 
pa radigm. I n so me wa ys our modeling is reminiscent of that 
by IOLA95I and IPP95L but we explore a much wider range 
of satellite orbit, mass and density, and therefore a greater 
breadth of context to test the character of possible resulting 
tidal disruption (§2). We return to simpler, more prosaic, 
one-component, MFL satellite structures because the previous 
studies have not completely explored the full parameter space 
possible within this paradigm and therefore (in our opinion) 
have not definitively ruled out MFL models. It is also possible 
that MFL dwarfs are a reasonable model for disrupting two- 
component models, after the o uter envelope of extended d ark 
matter has been stripped away (Klimentowski et al. 20071) . 

But a key differ ence between our work and the influential 
IOLA95I and IPP951 studies is that we explore the case of tidal 
disruption in dark matter dominated MFL systems. Also, 
in contrast to many previous efforts, our models are directly 



7 We note that Wilkinson et al. 1 2004) reported a sudden drop in projected 
velocity dispersion at about the King li miting radius for the Ursa Minor and 
Draco dSphs, and Kleyna et al. 1 2004) reported a similar finding for Sex- 
tans. Read et al. (2006b) attempted to model this apparent kinematical fea- 
ture and concluded that (1) dSph galaxies could have masses as high as 10 10 
Mq, and (2) tidal stripping is unimportant interior to ~ 1 kpc, although 
"tidal shocking" may st i ll be impo rtant. However, as has been shown by 
Lokas. Mamon & Prada (2005) and Munoz et al. 12005), the reported "cold 
points" are likely an artifact of binning and of ignoring "velocity distribu- 
tion outliers". In one case (Ursa M inor) the result was not confirmed with 
additional data ( Munoz et al. 200j). 
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compared to the most constraining, extant data set now avail- 
able, that for the Carina dSph. While Carina resembles most 
other dSphs both structurally and dynamically, it offers the 
advantage that it n ow has the most exte nsive spatial and kine- 
matical coverage (Munoz et al. 2006a) for any dSph outside 
the Sgr system. We specifically do not address the Sgr dSph 
because it has already been e xplored with similar methods and 
models bv lLaw et al.l d2005l) . and because Sgr has often been 
argued to be the exception, rat h er than the arch e type of Galac- 
tic dS phs (e.g., lMateolll998t iGilmord l2004t iGilmore et all 
2007). 



In IMufioz et al. (2006a) we presented empirical findings 
and arguments that Carina must be forming tidal tails, based 
on the appearance of the satellite at large radii (e.g., Fig. 1), 
the implied mass of the system if the outermost Carina stars 
are bound, and the observed velocity shear along Carina's ma- 
jor axis (Fig. 2d). The results of our modeling here lend con- 
siderable support to the hypothesis that the Carina system is 
presently tidally disrupting. We show that initially equilib- 
rium, Plummer model satellites on very shocked, radial orbits 
are transformed into systems with properties remarkably sim- 
ilar to Carina (and other dSph) satellites. In these models, the 
observed extended populations and the flat dispersion profiles 
are both a product of tidal debris, despite the presence of dark 
matter (in an MFL configuration). 

While we do not claim to have a unique model to describe 
the nature of the Carina system, our results dispel several 
common misconceptions about dSph galaxies and tidal dis- 
ruption often invoked to dismiss the tidal disruption hypoth- 
esis: We find that (1) despite the current trend of relying on 
ever more complex structural models (i.e. more free param- 
eters) to model dSphs, simple, one-component, tidally dis- 
rupting satellite models can adequately reproduce the strong 
physical constraints provided by existing observations of at 
least some dSph galaxies. (2) Our well-fitting tidal disrup- 
tion model for Carina even matches the central velocity dis- 
persion and M/L ~ 40 of the dSph. The presence of dark 
matter and tidal disruption are not mutually incompatible — 
MFL, dark matter dominated systems can produce tidal debris 
streams that are consistent with current observations. This 
has already been suggested b y our similar mo deling of Sgr 
dLaw et al.(l2005l) and Leo I dSohn et al.ll2007l) . but now we 
demonstrate the principle again with a high central M/L dSph 
system. Thus, attributing ever larger M/L and extended dark 
matter halos to dSphs to explain their observed outer structure 
and dynamics may not be necessary. (3) Despite considerable 
periodic shocking, significant tidal disruption, and the pres- 
ence of unbound stars near the dSph core, it is found the cen- 
tral velocity dispersion of the satellite still gives a reasonable 
estimation of its bound mass using the traditional core-fitting 
methodology. The influence of tides is primarily relegated to 
the oute r parts of the dSph, wher e the tidal force is strongest 
(see also Mash chenko et al.l2006l) . (4) The presence of even a 
small fraction of unbound stars can influence the observed dy- 
namics of dSphs at large radii. (5) Sgr may not be a unique ex- 
ample of a tidally disrupting dSph among present MW satel- 
lites. We conclude by explaining how tidally disrupting sys- 
tems that are well represented by an MFL paradigm are not 
necessarily at odds with the current A-CDM paradigm. 

An ancillary goal of the paper is to demonstrate a method- 
ology to navigate the parameter space of satellite orbit, mass 
and density to find matches to the observed structure and dy- 
namics of specific dSph galaxies having dynamical and struc- 
tural data to large radii using one-component, MFL, iV-body 
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FIG. 1. — Distribution of bona fide Carina RV-members on the sky from 
IMufioz et alJ2006al . 

simulations. In this case, we apply the methodology to the 
Carina system, but it is valid for application to any other sys- 
tem with sufficient structural and dynamical data. A second 
demonstration of the utilization of this methodology ha s been 
given in our study of the Leo I system dSohn et al.ll2007l) . 

2. W-BODY SIMULATIONS: MODEL PARAMETERIZATIONS 

Several properties of C arina we have o bserved 
dMaiewski et all 120001 120051: IMufioz et alj l2006al) — its 
substantially extended radial density profile, the highly elon- 
gated distribution of RV members, the flat velocity dispersion 
profile — are consistent with expectations for a satellite 
galaxy undergoing tidal disruption. But this interpretation is 
perhaps not readily acceptable or satisfying until a legitimate 
model or simulation can demonstrate that it is physically 
possible. Our goal here is to determine under what circum- 
stances (i.e. satellite mass, orbit and density distribution) 
— if any — the observed Carina velocity and structural 
characteristics can be imitated by simple, MFL satellites 
including tidal effects. We h ave modeled the Carina dwarf 
with 10 s particles initially in a Plu mmeridl911l) configuration 
and with the ensemble orbiting within a static, smooth MW 
potential. The self-gravitating satellite particles have their 
mutua l interactions calculated u sing a self-consistent field 
code (Hernquist & Ostrikerl [l992|). The form of the adopted 

a |]V 



MW potential, consisting of a (Miyamoto & Nagail (fl975l) 
disk, Hernquist spheroid and logarithmic halo, as well as the 
w orkings of the simu lation code is the same as that described 
in lLaw et all ((2005). The valu es for the adopte d potential 
parameters are those used in the lLaw et alj (|2005) simulation 
of satellite disruption that provided the 'best fit' to the spatial 
and velocity distribution observed for tidal debris associated 
with the Sagittarius dwarf spheroidal galaxy. In particular, the 
adopted halo potential is mildly oblate, with minor-to-major 
axis ratio of q = 0.9. 

To navigate the larger parameter space of satellite mass, 
density and orbit, we first ran a sampling "grid" of satel- 
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lite parameterizations to zero in on the family of models 
that best matches the observed Carina properties (§3.1). We 
then experimented with "fine tuning" models near the "best- 
matching" part of parameter space in an attempt to improve 
the match to Carina (§3.2). 

The properties of our initial baseline grid of models include 
the following considerations: 

(1) Current satellite position and velocity: The model Ca- 
rina is placed on orbits that yield a present Galactic location, 
heliocentric distance R and systemic RV for the satellite simi- 
lar to that observed for Carina (/ = 262°. 1, b = -22°. 2, /?= 101 
kpc and vqsr = 8.0 km s" 1 ). 

(2) Satellite mass: In our first tests Carina is represented 
with initial masses of either 10 6 , 10 7 , or 10 s M B . 

(3) Orbital plane: Based on previous (e.g.. lOLA95h and our 
own preliminary models of Carina-like systems that demon- 
strate satellite elongation by tides, we set the Carina or- 
bital plane by the position angle of the satellite (65° ± 
5°, or roughly 172° ± 5° i n Galactic coordinates, from 
llrwin & Ha tzidimitriou 1995). Whether or not this is true, 
Figure 1 shows that at least empirically the more extended 
Carina RV members do lie primarily along the central Ca- 
rina position angle. Preliminary simulations of satellites or- 
biting with a variety of orbital poles were run, and two (op- 
posite) orbital poles were found to satisfy the generation of 
debris at the observed PA, namely models with a pole of 
(l,b) = (159.3,-12.9)° or its antipode at (339.3, 12.9)°. Given 
the symmetry of our model potential and the likelihood that 
Carina is at its apoGalacticon (see below), we expect similar 
results with model satellites run on orbits with either of these 
poles, and elect to run the bulk of our model satellites on or- 
bits with the first pole. Once we find a favored model with 
this pole, we run a model with the same overall parameters 
but with the opposite pole (§3.2) to verify the symmetry of 
solutions. At this poi nt we do not invok e the measured proper 
motion of Carina by Piat ek et all ((2003), since the relative er- 
rors in this measurement are large and do not well constrain 
the direction of motion of Carina (the 95% confindence inter- 
val allows a nearly 80° spread in orbital inclination); however 
we point out some general consistencies of our best fitting 
model with Piatek et al. results in §3.2.3. 

(4) Orbital shape: Because the Carina velocity in the Galac- 
tic Standard of Rest convention (vgsr) is nearly zero, Carina 
must presently be near apogalacticon, near perigalacticon, or 
be on a nearly circular orbit. Given that the satellite radial 
velocity must be fixed to that observed, the ellipticity of the 
model orbits is set by the adopted transverse velocity, v, mns 
for a satellite at the given distance. We have created models 
characterized by v par = v, mns /v C i rc from 0.3 to 1.0, where v c ; rc 
is the Galactic ci rcular velocity at the Solar Circle, taken to be 
210 km s" 1 (see lLaw et al.ll2005h . The parameter v e ,, r sets the 
scale of the Galactic mass profile according to the adopted po- 
tential. The resulting orbits correspond to (peri/apogalactica) 
of (15/103) kpc to (85/103) kpc for orbits where Carina is 
presently at apogalacticon (orbits 1-7 with v par < 0.8). Or- 
bit 8 with v par = 0.9 is relatively circular and orbit 9 with 
v par =1.0 has Carina near perigalacticon (in an orbit with a 
152 kpc apogalacticon). Table 1 summarizes the main proper- 
ties for these orbits, including peri- (R per i) and apogalacticon 
(Rapo) distance, radial period (T mc i), range of peak velocities 
observed at perigalacticon (v pm ), the timescale of the tidal en- 
counter at perigalacticon (T enc = R pe ri/v per i), the mass of the 
MW interior to perigalacticon radius (M per i) and the Galactic 
density (p per d at that radius in our potential. We attempted 



orbits with v par as low as 0.2 ("orbit 1 "), but these have very 
small perigalactica that cause the orbit to be scattered by the 
disk potential so that no "typical" orbital properties can be 
given. 

(5) Simulation length and time sampling: For the initial set 
of models, Carina's orbit is traced backwards with a single 
point mass for the lesser of either five radial orbits or 10 Gyr, 
and then is evolved forward from this point in full A^-body 
mode. We have run most models with 4000 even time steps, 
but have found that the high perigalacticon velocities for some 
of the more radial orbits led to conservation of energy prob- 
lems at perigalacticon, and finer time sampling was imple- 
mented as necessary. Energy was conserved at the ~ 1 % level 
(of the internal satellite energy) when orbit 2 models were run 
with 64000 time steps, when orbit 3 models were run with 
16000 time steps and when all rounder orbits were run with 
4000 time steps. 

(6) Satellite density: Apart from the mass, we vary the 
physical scale length, ro, of the m odel satellite d istribution, 
which is initially set to generate a lPlummerl (1191 ll) model: 

$ = _ GM Carina ,Q (1) 

\A" 2+r o 

The density of the satellite, critical to the degree of robustness 
of the satellite against Galactic tidal forces and therefore a pri- 
mary driver in the satellite disruption rate, is set by the initial 
mass, Mcarina.o, ar, d scale, ro of the satellite. To maintain a 
constant mass density, the physical scale is sized nominally 
as 

r = 0.9(M Canna ,o/10 9 M ) 1 / 3 . (2) 

However, to vary the disruption rate, we adjust the scale by 
up to a factor of five (i.e., density varied by more than a factor 
of 100) at each adopted satellite mass. In all cases, the initial 
Plummer-configured satellite is allowed to equilibrate outside 
the MW potential before being introduced to the potential at 
the apogalacticon of the orbit. 

In summary, our initial grid of 80 models has been con- 
structed with three different satellite masses, with three or 
four different scale lengths per mass, and with eight differ- 
ent orbital shapes (orbits 2-9). Further explorations of satel- 
lite parameter space (§3.2) include additional combinations of 
mass, scale, number of radial periods run in the potential, and 
number of time steps in the integration. These extra models 
pertain primarily to orbit 2 and 3 models (which are found 
to give the best matches to the data, as shown below). In 
the end, nearly 200 different A^-body simulations to model the 
Carina dwarf galaxy were run on a nearly devoted Sun Blade 
1500 workstation, requiring approximately 6, 24, or 96 hours 
to run simulations having 4000, 16000 and 64000 time steps, 
respectively. Table 2 summarizes most of the combinations 
we have run by scale and mass, with each table entry indi- 
cating which orbit shapes were run for a particular mass and 
scale (and with the bold-faced entries representing the initial 
parametrization grid of 80 models). Asterisks indicate when 
a combination of mass, scale and orbit were run for different 
numbers (e.g., 1, 2, 3, 4, 5) of total radial orbits. 

3. RESULTS 

3.1. Results from the Initial Grid of Satellite 
Parameterizations 

For analysis of our initial grid of simulation parameteriza- 
tions we adopt as primary discriminants for the "success" of 
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TABLE 1 
Properties of Model Orbits 



orbit Vpar Rperi (kpc) R apo (kpc) T md (Gyr) v^fkms ') T mc (Myr) M pf ,„ (M Q ) p,,„, (M kpc 3 ) 



1 


0.2 


< 10 


> 145 


(not stable) 


2 


0.3 


15 


103 


1.46 


3 


0.4 


22 


103 


1.56 


4 


0.5 


34 


103 


1.67 


5 


0.6 


48 


103 


1.85 


6 


0.7 


64 


103 


2.05 


7 


0.8 


85 


103 


2.28 


8 


0.9 


95 


124 


2.71 


9 


1.0 


103 


152 


3.06 



a model its ability to match simultaneously (1) the observed 
velocity dispersion trend as a function of major axis and ra- 
dial distance from satellite center, (2) the general shape and 
size of the Carina radial density profile, and (3) the mean pro- 
jected velocity along the major axis of Carina. Fig ure 2 shows 
these o bserved constraints, with data taken from Mun oz et alj 
(l200oah . 

However, before moving directly to the best fitting mod- 
els for Carina (see §3.1.4), we take advantage of the fact that 
our systematic exploration of satellite parameter space leads 
to results that are useful in a more general context, revealing 
trends in the nature of satellite tidal disruption that are helpful 
for understanding the observational properties of tidal debris 
and disrupting satellites. To keep the scope of this more gen- 
eral discussion to a managable level, we limit our discussion 
primarily to the nature of the models as they appear at the dis- 
tance of Carina, which is at apogalacticon for orbits 1-8 and 
at perigalacticon for orbit 9. We make this simplification be- 
cause Carina's distance is more or less typical for the classical 
Galactic dSphs, and satellites on elliptical orbits spend more 
time near apogalacticon than in any other orbital phase. 

3.1.1. General trends 

To illustrate the general properties of the model satellites in 
our initial grid we have elected to show, first, the salient fea- 
tures of three satellites of the same mass, 1 x 10 6 M©, but of 
different scale lengths (or, equivalently, densities) evolved in 
our Galactic potential. As we have done for the real Carina 
data in Figure 2, we show in Figure 3 for the models with two 
different densities (a) the projected number density profile, 
the velocity dispersion profiles as a function of both (b) radial 
distance from the center of the satellite and (c) position along 
the major axis (normalized to r/,„,), which we characterize as 
a function of Galactic latitude, b (since this is a good approx- 
imation to the major axis for Carina's position angle), and 
(d) the mean projected velocity along the major axis. These 
properties are shown for four different orbital shapes, from 
very eccentric (v par = 0.3) to rather circular (v par = 1.0). In 
all models at least some degree of tidal disruption occurs and 
tidal debris dominates the structure and dynamics at large ra- 
dius. 

Some primary sweeping trends are demonstrated in Figure 
3. For example, the central velocity dispersion remains al- 
most unaltered by tides (in the sense that its value does not in- 
flate) until the satellite is near complete destruction whereas, 
at large radii, the velocity dispersion highly correlates to the 
de gree of ti dal stripping (see also Fig. 6). In fact, as found 
by IOLA95I, as long as the satellite retains a predominantly 
bound core, the value of the central velocity dispersion corre- 
lates with the instantaneous bound mass of the satellites. In 
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32.6 


1.90x10" 


2.47 x10 s 


310-380 


55.3 


2.41x10" 


1.22x10 s 
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92.8 


3.29x10" 
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1.32x10 s 
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357.8 


6.81x10" 
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198-222 


408.9 


7.48x10" 


5.91 xlO 4 


218-236 


417.0 


8.02x10" 


5.01 xlO 4 



the cases where the system is completely disrupted (i.e., indi- 
cated in the figure by having 0% bound mass today), a totally 
flat velocity dispersion profile is observed, and the magnitude 
of the dispersion is no longer correlated to the starting mass, 
but is also influenced by the initial satellite density and orbital 
shape. Density distributions, on the other hand, seem to be, 
for the most part, only mildly affected by the degree of tidal 
disruption when compared to the velocity dispersion profiles 
for a given mass/density/orbit combination. The main effect 
of tidal disruption on the density profiles is for the dSphs to 
develop an extended component that "breaks" from the ini- 
tial Plummer config uration, as has already been found in pre- 
vious studies (e.g.. [Johnston etal.lll999bi iRead et al.| [2006b; 
Klimentowski et al. 1 120071) . An overall enlargement of the 
satellites occurs only when the systems are near destruction, 
after which they lack a core in the density profile and present 
a rather flat spatial distribution. 

The general appearance of the debris for low eccentricity 
orbits is illustrated in Figure 4, where we show the Xqc-Yqc 
and Z gc -Yq C Cartesian projections (centered on the Galactic 
center) for a model with an initial mass of 1 x 10 7 M Q in two 
different, rather benign orbits. Tidal debris released at differ- 
ent times is represented by different colors, with the yellow 
particles marking the presently bound population. The upper 
and lowel panels show the model in a v par = 0.8 (nearly cir- 
cular) and = 0.5 (3:1) orbit respectively. In both cases, 
rather narrow, coherent debris streams are formed. 

In constrast, Figure 5 shows the Xqc - Yqq an d Ycc ~ Zee 
Cartesian projection of two 1.0 x 10 7 M Q satellites of differ- 
ent densities in a very radial orbit (v par = 0.3 and R per i = 15 
kpc) after five perigalacticon passages. The upper panels 
show a model with a scale of ro = 0.280 kpc (model 142) 
and the lower ones a model with ro = 0.194 kpc (model 102) 
— i.e. lower and higher density respectively. In both cases 
the debris is more broadly dispersed, but while model 102 re- 
tains a predominantly bound core after five orbits, model 142 
has been completely disrupted and lacks a core in the density 
profile (a property observed in all of our model satellites that 
have been completely destroyed). As Figures 4 and 5 demon- 
strate, for a given orbit and mas s, the mass-loss-rate is driven 
by th e density of the satellite (Johnston, Hernquist, & Boltd 
Il996t see also Figure 6, where the density is kept constant 
and the mass-loss-rate stays the same for a given orbit, even 
as mass is varied). 

Other basic trends depend heavily on the orbit shape, e.g., 
more circular, rather gentle orbits (apo/perigalacticon ratio no 
greater than 3:1) versus more energetic, rather eccentric or- 
bits. These results are summarized in Figure 6, which shows 
the density, mean velocity and velocity dispersion trends for 
three disrupting satellites of different masses (1 x 10 6 , 1 x 10 7 
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FIG. 2. — (a) Number density profile for the Carina dSph. (b) Observed 
velocity dispersion as a function of radial distance from the Carina's center, 
normalized to /-/„„ .(c) Observed velocity dispersion as a function of Galactic 
latitude, b, which is a good proxy for the major axis of Carina. We shift the 
abscissa to be centered on the Galactic latitude of the center of Carina, b c . (d) 
Observed mean projected velo city as a function of major axis distance. All 
data from Munoz et al. (2006a). 



and 1 x 10 8 M Q ) in two different orbits (v par = 0.8, left pan- 
els; v par = 0.3, right panels), this time with central density kept 
constant (by adjusting the satellite scale to track the cube root 
of the mass according to Eq. 2) in order to explore the trends 
with mass. 

After five orbits (or, alternatively, 10 Gyrs, see §2.1), most 
of the initial mass still remains bound for model satellites in 
gentle orbits (left panels of Figure 6), regardless of their ini- 
tial masses. At the same time, all of the simulated satellites 



in these orbits experienced some degree of tidal stripping that 
leads to extended spatial components populated by tidal de- 
bris. The left panels in Figure 6 show that model satellites 
with large perigalacticon distances and almost circular orbits 
(R P eri=85 kpc and v par = 0.8) undergo relatively modest tidal 
stripping and their velocity dispersion trends are found to de- 
cline over degree (kpc) scales, as expected for mostly bound 
systems (containing particles on predominantly non-circular 
internal orbits) where the velocity dispersion is the highest at 
the center and then decreases at large projected radii. How- 
ever, the velocity dispersion at large radii still departs from 
the initial, unperturbed value, indicating some tidal heating of 
the bound population at large radii. 

As the orbits become less circular (and the perigalactica 
smaller), tidal stripping becomes more vigorous and more un- 
bound particles are produced. The right panels in Figure 6 
show how tidal debris (unbound particles) generally produce 
flat or rising dispersion profiles at large radii from the dSph 
center. In fact, the size of the RV dispersion in the tails corre- 
lates with both the mass of the satellite and the orbital shape as 
long as the satellite retains a predominantly boun d core; this 
is as w as found in the modeling of the Sgr system dLaw et al.l 
2005), and was, in fact, one way that was used to constrain 
the mass of the Sgr system in that study. 

A key feature demonstrated by these radial orbit models — 
indeed by all models retaining bound cores (e.g., in Figures 3 
and 6) — is that the central velocity dispersion is still domi- 
nated by the bound population even when the outer parts have 
greatly inflated velocity dispersions. Thus Galactic tides have 
little effect on the inferred central M/L, even in systems with 
high disruption rates. However, a radial rise to a higher veloc- 
ity dispersion is also a characteristic product of the unbound 
component in these radial orbit models. This rise in veloc- 
ity dispersion can happen well within the limiting radius of 
the bound system, reflecting the presence of unbound parti- 
cles quite near the center of the satellite. To illustrate this, 
in the left panels of Figure 7 we show the spatial distribu- 
tion (in Galactic coordinates) of particles near the limiting ra- 
dius for a model in an eccentric orbit (1 x 10 7 Mq, ro = 0.194 
kpc, Vp ar = 0.3). Panels from top to bottom show particles 
that have become unbound in successive orbits, with the top 
panel showing the presently bound population. This figure 
demostrates that, in general, the centralized tidal debris re- 
sponsible for inflating the velocity dispersion radially, is dom- 
inated by the "puff" of debris released on the last perigalactic 
impulse, although particles unbound from prior perigalactica 
can linger near the dSph center for many orbits. For compari- 
son, the right panels of Figure 7 show the same feature but for 
the model satellite in a more benign orbit (v par = 0.8) where 
the tidal debris dissipates much more quickly. 

3.1.2. Tidally Induced Rotation? 

PP95 pointed out that the most unique and readily observ- 
able consequence of tidal disruption is tidal shear that resem- 
bles solid-body rotation along the major axis of the dSph. To 
compare our results to these expectations, the lower panels of 
Figures 3 and 6 show the mean projected velocity along the 
major axis distance for the same models shown in the upper 
panels. Satellites in nearly circular orbits show, in general, 
a smooth velocity gradient within their r/ lm (e.g., Fig. 6, left 
panels). However, this small velocity gradient is not a result 
of Galactic tides (there is only a small degree of tidal disrup- 
tion suffered by satellites in these orbits) but is predominantly 
a result of the line-of-sight projection of the satellite mean ve- 
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FIG. 3. — From top to bottom: Projected number density profiles, velocity dispersion profiles as a function of radial distance and distance along the major axis 
(characterized as a function of Galactic latitude and normalized to the respective r/„„) and mean projected velocity along the major axis for {left panels) model 
satellite of mass 1 X 10 6 Mq and scale length of i'o=0.090 kpc, in different orbits (v par = v tmns /v c i rc =0.3, red lines; 0.5, green lines; 0.7, yellow lines and 1.0, blue 
lines), and (right panels) model satellite of mass 1 X 10 6 Mq in the same orbits as the left panels, but this time with scale length of ro=0. 194 kpc. The percentage 
numbers represent the fraction of bound mass remaining at the end of the simulation. Solid lines represent all particles while dotted lines show the behavior of 
the bound particles only. 



locity along the orbital path. On the other hand, satellites in 
m ore r adial orbits do more strongly show the effect described 
bv lPP95l (e.g., Fig. 6, right panels). While satellites that have 
been completely destroyed show, for the most part, a smooth 
gradient along the major axis, satellites on more radial or- 
bits that retain at least some bound core show a "rotation" (S- 
shape) signature like that described by PP95. However, it is 
important to note that this signature is observable mostly out- 
side the r/,„, of the satellite core, whereas within r/„„ the mag- 
nitude of the signal is typically relatively modest — e.g., small 
compared to the central velocity dispersion. In these cases of 
satellites with bound cores, the bound component shows al- 
most no signature of tidal shear and the "rotation" is reflected 
only in the unbound particles (thus, this "shearing" or "rota- 
tion" is really just a reflection of how the tidal debris with its 
net angular momentum and energy is organizing in reponse to 
the MW potential). These results suggest that, while apparent 
rot ation c an be a consequence of tidal stripping, as suggested 
by IPP95L for satellites that are not completely destroyed this 
effect can occur well outside the typical radii to which most 



current R V surv eys of dSphs have been limited. Thus, in most 
cases the IPP95I signal cannot be counted on as evidence for 
lack of tidal effects since the rotation would generally not be 
expected to be seen, or, if detected, it would likely be deemed 
as "statistically insignificant" in most studies. We return to 
this point in §4.1. 

3.1.3. Time Dependence of the Satellite Properties 

As mentioned in §1.1, a primary problem with previous 
modeling of tidally disrupting, low M/L satellites is that any 
successful matches to observed dSph satellites depends on 
catching these systems in one particular phase of evolution. 
However, in all cases we have explored, save those where the 
satellite is completely destroyed, the radial density profile of 
the satellite resembles an "inner King+extended population" 
profile that looks nearly the same for different orbital passages 
when viewed at the same orbital phase (uppermost panels of 
Figures 2, 3, 6, 8 and 9). Figure 8 demonstrates our general 
finding that the overall spatial and dynamical character of dis- 
rupting satellites can be maintained for long periods of time 
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FIG. 4. — Xgc~Ygc and Ygc~Zgc Cartesian projections for a model with 
a starting mass of 10 7 Mq in two rather benign orbits. The upper panels show 
the satellite in a v par = 0.8 orbit while the lower panels show the system in a 
Vpar = 0.5 orbit. The coloring scheme is as follows: Yellow dots mark parti- 
cles that are still bound, while green, blue, red, cyan and magenta represent 
material that became unbound after 1,2,3,4 and 5 orbits respectively. 



— as long as the satellite maintains a bound core but under- 
goes a more or less steady mass-loss-rate. This is an advan- 
tageous property of our models over the previous, more spe- 
cialized A^-body simulations of tidal d isruption. In contrast 
to, for example, the findings of IPP95L the specific physical 
satellite properties we seek in the models are not rare phe- 
nomena that we have to "catch" at special times in the life of 
the satellite, but rather long-lived phenomena that require no 
invocation of particular timing (apart from observing the sys- 
tem at the proper orbital phase — near apogalacticon in the 
case of Carina, which is a position near where satellites spend 
a considerable fraction of a radial perio d). Our dSph mod- 
els also differ from those of lKuhnl (1 19931) and lKroupal (1 19971) 
where the satellites must be observed at very specific times 
in their tidal evolution (namely, after total disruption or just 
prior to total disruption for these two references, respectively) 
and along specific lines of sight. 

On the other hand, Figure 8 shows that while the overall 
character of the satellite density and velocity dispersion pro- 
files is maintained for long periods of times, there is some 
gradual evolution in the scaling of the density profile and the 
magnitude of the velocity dispersion at large radii that is de- 
pendent on the amount of past disruption. Figure 8 shows the 
velocity dispersion profile observed today for satellites with 
the same initial mass (in this case 10 6 Mq) on the same radial 
orbit, but which have been dropped into the Galactic poten- 
tial at successively earlier apogalactica. In this case, a central 
dip to a progressively lower central velocity dispersion is ob- 
served, while the dispersion of the debris-dominated part of 
the system changes with time, with slighly colder debris ac- 
cumulating after the satellite has suffered more mass loss. Of 
course, this mass loss causes the overall density of the satellite 



FIG. 5. — Xqc-Ygc and Ygc~Z<x Cartesian projections for two mod- 
els in rather eccentric orbits (v par = 0.3) after five perigalacticon passages. 
Both models had an initial mass of 1 X 10 7 Mq but the upper panels show 
the model with a rg = 0.28 kpc (less dense) while the lower panels show the 
model with a ro = 0.194 kpc (more dense). The coloring scheme as the same 
as in Figure 5. As may be seen, the lower density model is destroyed af- 
ter five orbits (no yellow component) but the higher density model keeps a 
predominantly bound core. 

(at all radii) to drop with time. The trend of the phenomena 
shown in Figure 8 thus admit yet another "free parameter" — 
i.e. number of orbits of tidal interaction — that one may tune, 
in addition to orbit shape, initial satellite mass, and satellite 
scale, to match the properties of any particular real satellite, 
such as the Carina system. 

3.1.4. What is needed to match Carina? 

In only one part of the systemic survey of orbit, mass, den- 
sity parameter space do the simulations reproduce the general 
character of the observed Carina velocity dispersion profile as 
well as the general projected density distribution of Carina. 
The operative conditions in common among the "successful" 
models is a satellite that experiences intense, periodic tidal 
shocks at the perigalacticon of a very radial orbit, but one 
that also retains a bound core. Most of the models run with a 
Vpar = 0.3 satellite orbit (orbit 2), which have perkapogalactica 
of 15:103 kpc and radial periods of 1.5 Gyr, share this char- 
acter (e.g., Fig. 8), as do some model satellites on the less 
extreme orbit 3 (v par = 0.4). 

We may summarize the results of this initial survey of pa- 
rameter space to find a model that resembles Carina as fol- 
lows: 

• It is possible to find MFL iV-body models of disrupt- 
ing satellites, constrained to the present position and 
radial velocity of Carina, that reproduce the general 
morphological character (i.e., a "inner King+extended 
population" radial density profile) and dynamical char- 
acter (i.e., a flat/rising dispersion profile with radius) 
observed in the Carina system. 
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FIG. 6. — Same as Figure 3 but this time, (left panels) model satellites of masses 1 X 10 6 (red lines), 1 X 10 7 ( green lines) and 1 X 10 8 Mq (blue lines) with 
initial central density kept constant by adjusting the Plummer scale length (see text) in a rather circular (v par = 0.8) orbit, and (right panels) the same model 
satellites but now in a eccentric orbit (v par = 0.3). As in Figure 3, the percentage numbers represent the fraction of bound mass remaining at the end of the 
simulation. Solid lines represent all particles while dotted lines show the behavior of the bound particles only. 



• The common properties of models succeeding in this 
way are that the satellite retains a bound central core 
but experiences periodic, impulsive shocks on a rela- 
tively radial (e.g., 7:1) orbit that spawns nearby clouds 
of unbound tidal debris. Carina resembles the system as 
observed near the apogalacticon of one of these orbits. 

• The rise in the velocity dispersion away from the core 
is caused by the growing contribution of unbound de- 
bris which becomes important even well inside the pro- 
jected limiting radius of the bound component. The ris- 
ing dispersion profile at larger radii is due to the domi- 
nance of tidal debris there. 

• Our simulations do not depend on special timing to 
match the general Carina properties, other than observ- 
ing the system near the correct orbital phase. Thus, 
the results of our modeling do not hinge on particularly 
special circumstances other than that the satellites be on 
rather radial orbits. 

In the next section we explore this family of models further, 
both to find a better match to the Carina data by fine-tuning 



the orbital properties, and to understand the implications of 
such models for the possible physical state (e.g., mass, M/L, 
mass-loss-rate) of the Carina system. 

3.2. Fine Tuning the Models 

3.2.1. Zeroing in on Carina 

>From the previous general survey of single component N- 
body simulations we understand that orbital shape plays a crit- 
ical role in producing the dynamical effects we seek to emu- 
late in the Carina system, and that very radial orbits are pre- 
ferred for a good match to the observed Carina properties. To 
obtain an even better match to the observations, it is worth- 
while investigating a finer gradation of parameter space for 
these radial orbit models. For now, we focus on models with 
satellites on orbit 2. With our single-component models this 
leaves as free parameters to vary: the initial satellite total 
mass, scale, and duration of the numerical integration (i.e., 
number of radial orbits). These parameters are expected to 
affect the satellite and debris in certain ways for a particular 
orbit: 

• The measured King core radius (r core ) and r/„„ will ob- 
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TABLE 2 
Model Parameterizations" 



Scale 



initial mass (Mq) 


0.045 


0.090 


0.194 


0.280 


0.350 


0.418 


0.550 


10 6 


2 


(l-6)*,7-9 


2-9 






2*,3-9 




3xl0 6 


1,2* 


1,2* 












10 7 




(l-4)*,5-9 


2*,3-9 


2*,3-9 




2-9 




3 X 10 7 








2 








3.75X10 7 










2 






4.5 xlO 7 












2 




10 8 




2-9 


2-9 






1-9 


2 



Entries give orbits run for each satellite and scale mass. Boldface entries highlight the initial parameterization grid. Asterisks denote models run with multiple versions having the 
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FIG. 7. — (a) Spatial distribution of bound particles for a model of 1 x 10 7 
Mq and ro = 0.194 kpc. Panels (b), (c) and (d) show the spatial distribu- 
tion of particles lost during the last and successively previous orbits respec- 
tively (equivalent to the different colors debris shown in the other models 
portrayed in Figures 4 and 5). Left panels show the satellite in an eccen- 
tric (v por =0.3) orbit while the right panels show the same satellite in rather 
circular (v par =Q.&) one. As may be seen, the nearby debris is dominated by 
particles released during the last orbit, although debris released in previous 
passages can linger near the center for many orbits in the case of eccentric 
orbits. 

viously be influenced by the input physical scale of the 
satellite, r . 

• The central velocity dispersion of the satellite is ex- 
pected to vary dimensionally as (mass/scale). 

• The amount of debris is set by the mass-loss-rate and 
this is correlated to the density of the satellite (for a 



given orbit). Thus it scales as (mass/scale 3 ). This will 
affect not only the relative density of the extended pop- 
ulation to the central density of the satellite in the radial 
profile, but, of course, how quickly the velocity disper- 
sion switches from being dominated by bound versus 
unbound debris. 

• As we hav e seen above and in our previous studies of 
Sgr debris (lLaw et al.ll2005h . the velocity dispersion of 
tidal debris correlates with the satellite mass as long as 
the satellite maintains a bound core. 

Guided by these rules of thumb, we attempt to vary the mass 
and scale of the model satellite in a way to best match the 
radial velocity dispersion and radial density profiles simulta- 
neously. 

First, we set the physical scale. Since we are looking for a 
model that matches all of the Carina properties, we start by 
requiring the model satellite to have a physical scale match- 
ing as closely as possible that of Carina. Considering that the 
size of the bound component of the model satellites does not 
vary significantly over their lifetime (e.g., Figure 8), we set 
the Plummer physical scale of our model Car ina guided b y 
the observed Carina core radius, r core = 258 pc (lMateolll998h . 
Not surprisingly, King distributions fitted to satellites mod- 
eled in our initial general survey (see Appendix A) show that 
those with ro = 0.280 kpc yield the closest match to the Carina 
density profile — i.e., those for which ro is similar to Carina's 
r 8 

1 core- 

Once the physical scale is set, we seek a model that is capa- 
ble of reproducing the observed c entral velocity dispe rsion of 
Carina (6.97 ± 0.65 km s _1 , from lMufioz et alJl2006ah . Since 
this varies as (mass/scale) and the scale of the model is already 
set, that leaves us with the initial mass as the pertinent vari- 
able. But by varying the initial mass we are also simultane- 
ously changing the initial density of the model and that drives 
the mass loss rate, which, in turn, affects the shape of the 
velocity dispersion profile as well as the current bound mass 
of the satellite; these observational constraints each tightly 
and independently constrain the choice of mass. Nevertheless, 
a model reasonably matching all of these constraints can be 
found. 9 For example, Figure 9 shows the velocity dispersion, 

8 While r CO re radius and the Plummer length scale of a given density dis- 
tribution are not strictly equivalent, their ratio is of order unity for low con- 
centration systems like dSphs. 

9 The fact that all of these constraints can be simultaneously satisfied lends 
considerable weight to the viability of the single component models; after all, 
were they an incorrect description of the physical state of Carina, it is easy 
to imagine our models having difficulty satisfying all observable constraints 
with the limited number of free parameters. 
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FIG. 8. — The evolution of the density, velocity dispersion and mean veloc- 
ity profiles over five orbits for one of the low mass, medium density, radial 
orbit models showing the general properties observed in Carina. This shows 
that satellites can retain the general character of a rising velocity dispersion 
over many orbits, but the asymptotic dispersion value is a function of the 
amount of past disruption. 

projected mean velocity and density profiles for satellites in 
a v par = 0.3 orbit with initial masses of 2 x 10 7 , 3 x 10 7 and 
1 x 10 s Mq, after 5 orbits. The lower and higher mass models 
yield current central velocity dispersions that are too low and 
high compared to Carina respectively, while the intermediate 
mass model gives the closest match. This intermediate mass 
model provides a better-matching density distribution as well. 

While Figure 9 demonstrates that of the three models shown 
the 3 x 10 7 M Q model provides the closest match to the shape 
of the velocity dispersion profile, its values are still systemat- 



FlG. 9. — Three models with an initial physical scale of ro = 0.280 kpc, 
chosen to yield present satellite sizes comparable to that of the luminous com- 
ponent of Carina (see text for details). Three different initial masses of 2,3 
and 10 x 10 7 Mq are shown to explore the resultin g number density, ve loc- 
ity dispersion and mean velocity profiles. Data from Munoz et al. 1 2006a) are 
shown for comparison as open symbols with error bars. 



ically low compared with the data. Further fine tuning of the 
model Carina thus requires an overall increase of the veloc- 
ity dispersion. We achieve this by increasing the initial mass 
of the model while at the same time making the orbit slightly 
more eccentric to keep a similar mass loss rate. The net ef- 
fect is a model that yields the right central velocity dispersion 
while at the same time providing a good match to the shape of 
the velocity dispersion and the number density profiles (those 
shown in Figure 10). 
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Finally, we note that models initially run and presented 
so far were designed to have an orbital pole of (l,b) = 
(159.3,-12.9)°. But when we eventua lly collected spectro - 
scopic data at large radius (presented in Mu noz et al.ll2006a) . 
we found that these models were actually yielding the op- 
posite mean velocity trend from that being found in Carina 
(Fig. 2). To fix this difference with our completed simula- 
tions, in the end we re-computed our best matching models 
by placing the satellites on orbits with the polar antipode of 
(l,b) = (339.3, 12.9)° to obtain the correct velocity trend. 

Tables 3 and 4 give the parameters for two models (502G 
and 5021) that yield close matches to the data. As may be 
seen, these "best matching" models to the Carina structure 
and dynamics have initial masses of 3.6 and 3.8 x 10 7 M , but 
have evolved after five orbits to present masses of ~ 1 .9 x 10 7 
Mq. Interestingly, these final masses are comparable to the 
current mass of Carina derived using th e core-fitting method 
dMateo et all 1 19931: iMufioz et al J[2(K>6:i: see §3.2.2). Figure 
10 compares the final state of these simulations against the 
data for Carina. In the top panels, we compare the number 
density profiles. In the middle pairs of panels we compare 
the velocity dispersion profiles both as a function of angular 
distance from the center of Carina and as a function of Galac- 
tic latitude respectively, and in the bottom panels we compare 
the mean velocity trend of the satellite in Galactic Standard 
of Rest velocities. For both simulations 502G and 5021, the 
models lie within 1-cr of the data for almost all positions in 
the dSph for all four of the observational trends shown. 

3.2.2. Mass-to-light ratios 

Tables 3 and 4 also give for the different models the evolu- 
tion of the satellite from its initial conditions at each apogalac- 
ticon and at the final state, with the distance, best fitting 
central King profile parameters (see Appendix A), fractional 
bound mass remaining, and central velocity dispersion given. 
As may be seen, the actual King profile radii of the satellites 
remain relatively constant through time. However, the central 
velocity dispersion is seen to drop as the satellite loses bound 
mass, closely tracking the instantaneous mass — as expected 
if the steady state virial theorem remains a good approxima- 
tion in the face of relatively modest mass loss rates. 

If we desire the M/L of the satellite to remain constant in 
time (i.e., dark matter and luminous matter are lost in equal 
proportion for MFL models), then the central surface bright- 
ness of the satellite must evolve scaled to a 1 . To demonstrate 
that this actually happens, we also give for each apogalac- 
ticon time stamp in Tables 3 and 4 a measurement of the 
global M/L as inferred using the standard core-fitting tech- 
nique (Illingworth 1976) applied with the velocity dispersion 
and core radius measured at that point in the simulation: 



(M/L) tot = 



166. 5/^/1 



(3) 



Application of this equation requires calibration of the models 
by the actual central surface brightness in order to "paint" an 
appropriate fraction of the particles as luminous. For the fi- 
nal time stamp of the models, which are intended to represent 
the present Carina dSph, we simply ado pt the observe d lumi- 
nosity of Carina of 4.3 x 10 5 L Q from iMateol dl998l) . Once 
this is done, we can determine the scaling from particles to 
luminosity by the ratio of the actual luminosity to the total 
number of particles in the final time step of the model. This 
scaling allows us to track the measured M/L of the satellite 



as determined through core-fitting backwards over time (e.g., 
the ninth column in Tables 3 and 4). But because we are as- 
suming MFL, we can of course independently assess what the 
actual M/L of the bound system is (which is constant by de- 
sign) by the fact that we know what the bound mass is at the 
last stage of the simulation and we are assuming the satellite 
has currently the luminosity of Carina. Tables 3 and 4 shows 
that the M/L derived by core-fitting reasonably matches the 
actual M/L of the model throughout its evolution. We have 
found that this condition breaks down only when the satellite 
reaches near c omplete dissolution, as previously reported by 
iKroupal (119971) . The stability of the observed M/L with time 
is another convenient feature of our models that frees us from 
concerns about special timing requirements to match models 
to data. Moreover, it is not a forced condition, but one that 
arises naturally in our models. 

3.2.3. Comparison to Proper Motion-Derived Orbit 

The orbital parameters of those models that give the best 
match to the observed Carina properties (Tables 3 and 4, Fig- 
ure 10), are in very good agreement with the one previous at- 
tempt to determine the orbit of Carina directly. Piate k et ail 
(2003) measured the proper motion of Carina using HST 
imaging and found that the dSph must currently be near 
apogalacticon, and they derive the best estimate for apo- 
and perigalactica at, respectively, 102 and 20 kpc; this is in 
very good agreem ent with our results of 103 and 15 kpc. 
Piatek et al] d2003l) also give an estimate for the orbital pe- 
riod of 1 .4 Gyr; again, this matches well our derived period 
of 1.5 Gyr. Although their study does not find Carina in a 
polar orbit (the inclination of their derived orbit with respect 
to the Galactic plane is 39°), the 95% confidence interval for 
their derived inclination is so wide (23°, 102°) that it includes 
a polar orbit such as we derived here 10 . 

In summary, our "fine tuning" procedures have allowed us 
to find MFL, tidal disrupting satellite models that quite suc- 
cessfully fit the prime observational properties of Carina: the 
velocity, velocity dispersion and radial density profiles. More- 
over, along with other properties that arise naturally, such 
as the close and correct tracking of the instantaneous central 
M/L, these model properties (1) do not require special condi- 
tions (e.g., particular viewing perspectives, critical evolution- 
ary states) to be observed, (2) have long term stability over 
the life of the satellite, and (3) do not require the invocation 
of extra free parameters beyond initial satellite mass, scale, 
and orbital shape. 

3.2.4. Tidally Induced Ellipticity? 

Despite being a good match for most of the structural and 
kinematical properties of Carina, the only shortfall of our 
models that we have foun d is that they fail to match Carina's 
observed ellipticity (0.33, Mateo 1998); rather, in general the 
models maintain bound components that are nearly spherical. 
Obviously, elongation of the model satellites is observed over 
very large scales, when tidal debris starts to organize into tidal 
tails, but not over radii where ellipticity is seen in the dSphs. 1 1 

10 Our own derived orbital inclination (§2) is based solely on the general 
orientation of the extended Carina population (Fig. 1 ; see also Fig. 3) and this 
estimate itself is likely subject to some 10° uncertainty. 

1 1 The orientation of the nascent tails in the simulations agrees with the 
observed instrinsic elongation of Carina, but of course the orbital poles for 
the simulations were chosen assuming that the major axis of Carina lies along 
its orbital path (§2). 
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FIG. 10. — From top to bottom: Number density profiles, velocity dispersion versus radial distance, velocity dispersion as a function of major axis, and mean 
RV trend (with velocities in the Galactic Standard of Rest fra me) for our two best matching models. Model 5021 is shown by the left panels and model 502G by 
the right ones. The data for Carina from Muiioz et al. 1 2006a) are shown as solid triangles. 



TABLE 3 
Properties of Model 502g 11 



orbit 


D (kpc) 


r CO re (pc) b 


r lim (pc) b 


Rtidal (PC)° 


cr„ (kms-') d 


M bomd (10 7 M ) 


Mbomd /M 


(M/L) cf ((M/L) a ) e 


(M/L),„ ((M/L) ) f 





104 


230 


725 


550 


8.73 


3.60 


1.00 


38 


45 


1 


96 


220 


710 


535 


8.48 


3.26 


0.91 


38 


45 


2 


97 


220 


760 


505 


7.96 


2.77 


0.77 


39 


45 


3 


106 


225 


770 


490 


7.60 


2.48 


0.69 


41 


45 


4 


97 


235 


700 


470 


7.22 


2.23 


0.62 


43 


45 


5 


101 


250 


815 


450 


6.83 


1.94 


0.54 


47 


45 


Carina 8 


101 


258 


846 




6.97 


1.76 




41 





While spherical cores are a general result for those models 
that retain a significant fraction (> 10%) of their initial mass 
as bound, satellites with a bound core but having the observed 
ellipticity and position angle of Carina are obtained only when 
~5% or less of the initial mass remains bound. Figure 1 1 
(upper panel) shows a coutour map of the density distribu- 
tion for a model that is near complete destruction. The elon- 
gated appearance of the satellite is clear. Interestingly, when 
they are near destruction, these satellit es exhib i t a kin ematical 
behavior similar to that described by Kroupi (1 19971) — that 



is, their central velocity dispersions are significantly inflated 
along with an associated artificial increase in the inferred cen- 
tral mass density and total mass. However, in these cases, the 
satellites no longer provide a good simultaneous match to the 
density and velocity dispersion profiles of Carina, as Figure 
12 shows, due, partly, to a sudden and significant inflation in 
the satellite size (by up to an order of magnitude) with respect 
to its starting scale length. We conclude, therefore, that if our 
MFL models provide the correct explanation for the structure 
and dynamics of Carina, the observed ellipticity in this satel- 
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TABLE 4 
Properties of Model 502i a 



orbit 


D (kpc) 


r CO re (Pc) b 


riim (pc) b 


Rtidal (PC) C 


o-vtkms- 1 ) 1 ' 


M bmmd (10 7 M ) 


M bou „ d /M 


(M/L) cf «M/L) ) e 


(M/L) m ((M/L) ) f 





105 


250 


900 


500 


9.00 


3.80 


1.00 


40 


43 


1 


97 


230 


740 


485 


8.70 


3.42 


0.90 


38 


43 


2 


97 


250 


770 


455 


7.89 


2.83 


0.74 


41 


43 


3 


106 


230 


715 


435 


7.59 


2.50 


0.66 


39 


43 


4 


97 


230 


760 


420 


7.16 


2.21 


0.58 


40 


43 


5 


101 


235 


850 


395 


6.60 


1.85 


0.49 


41 


43 


Carina 8 


101 


258 


846 




6.97 


1.76 




41 





a Run in an orbit with R a po/Rperi of 101/15 kpc. b Core and King radius from best fitting to bound population. Tidal radius from equation (2) ( Oh et al. 1992).^ Central projected 
velocity dispersion. e Global Mass to Light Ratio from Core-Fitting technique applied to the model's structural parameters/ Global Mass to Light Ratio from the known bound mass in 
the model. This value is constant over time by design. g D, r core and r^ m from Mateo i 1998). a v and (M/L) c f from Munoz et al. 2006a. 



lite is not likely tidally induced (see §4.2), and is probably 
related to the initial shape of the satellite. 

4. DISCUSSION 

Our work here has demonstrated several important features 
of tidal disruption and the dark matter content of dSphs in 
general, and of the Carina dSph in particular. These results 
are at odds with some commonly held perceptions regarding 
the observable effects of tidal disruption. These perce ptions 
hav e partly been driven by extrapolating the results of IPP95I 
and IOLA95I beyond the specific cases actually explored in 
those papers. We use our analysis of tidally disrupting TV-body 
satellites to readdress eight statements commonly made about 
dSphs and tidal disruption that we believe are either incorrect, 
or at least are overly simplistic. 

4.1. Misperception 1 : Lack of Rotation Signatures Means 
Tidal Disruption Is Not Happening 

Current RV surveys of Galactic dSphs extend now to the 
outskirts of these systems (as defined by the distribution of 
their luminous matter) and this provides new constraints use- 
ful for chec king th e dynamics of tidally disrupting TV-body 
simulations. IPP95I stated that the most useful result from their 
TV-body simulations is the identification of tidally induced ro- 
tation along the major axis of the system as a signature of 
tidal disruption. Subsequently, this has become perhaps one 
the most widely used criteria to assess whether tidal stripping 
is affecting particular dSph systems. 

To date, kinematical surveys of most dSphs have found lit- 
tle evidence for significant rotation. This lack of a clear ro- 
tational signature has thus lead to the conclusion that tidal 
strippi ng of Galactic dSphs is for the mos t part unimportant 
(e.g JKleyna et aLlbOOlHKoch et alJ2007allbf) . However, even 
careful examination of the PP95 results indicates that lack of 
rotation signatures within current RV surveys is not necessar- 
ily indicative of the importance (or lack thereof) of Galactic 
tides. IPP95I found significant rotation along the major axis 
of their model satellites only when the satellites were near 
perigalacticon of their eccentric orbits and when the host 
galaxy 's mass distribution was modeled as a point mass. But 
PP95 also noted that when the Galactic halo was modeled by 
a more realistic potential, this apparent rotation became less 
pronounced in their models and was only observed over large 
spatial scales (this point was also made by OL A95). Perhaps 
most relevant to the typical Galactic satellites, IPP95I saw that 
the rotation signal effectively disappears as the satellites move 
to apogalacticon, and this is an orbital phase near where satel- 
lites are most likely to be observed. Our more comprehen- 
sive exploration of satellite parameter space and modeling in a 



more realistic MW potential also indicates that tidally induced 
rotation is not always a good discriminator of tidal disruption. 
As discussed in §3.1.2, detectable levels of induced rotation 
along the major axis is only observed over scales larger than 
approximately r/,,„ of the model satellites, whereas most cur- 
rent RV surveys are confined to regions interior to these radii. 

Thus, because tidally induced rotation can be completely 
unobservable within the current spatial limits of most dSph 
RV surveys for satellites at typical points in their orbits, this 
test is not reliable for ruling out tidal disruption. Eventually, 
when large enough radii can be probed dynamically, this kind 
of test may be more appropiate; this is the case for Carina 
where RV members have been detected to at least 4.5r; lm and 
rotation along its major axis is observed beyond its ru m . 

However, an additional complication to using this "rota- 
tion" test over large angular scales is that even if a satellite 
is completely bound, eventually the varying projection of the 
satellite bulk motion on the radial velocity can produce an 
apparent "shearing" trend similar to the one caused by tidal 
effects. If the mean velocity as a function of position can be 
ascertained with great enough precision, it may be possible to 
discriminate the more uniform trend produced by projection 
from the more 5-shaped trend seen in many tidal disruption 
models, but the differences are subtle. 

4.2. Misperception 2: Tidally Disrupting Satellites Have 
Pronounced Ellipticities 

An interesting result from our general survey is that for 
models that retain a bound core with more than ~ 10% of 
the initial mass, the bound component of our satellites always 
remain spherical, with no perceptible, tidally induced elon- 
gation. Figure 1 1 (lower panel) demonstrates this situation: 
Here, we show the contours of the projected number density 
of one of our best matching Carina models only using parti- 
cles within two r/ im . Even though some debris is projected 
within r/„„, the overall azimuthal satellite density distribution 
is dominated by the still bound material and, as is readily ob- 
vious, the contours show no substructure or clumping and the 
satellite remains fairly round and regular. 

This lack of tidal elongation in our model seems to be, in 
principle, at odds with the findings of PP95 (see their Fig- 
ure 5). These authors found their simulated satellites enlarged 
and presented projected ellipticies as high as 0.7, which PP95 
ascribed solely to the effects of Galactic tides. But a care- 
ful comparison between our simulations and PP95's indicates 
that there is no actual disagreement between our models and 
theirs, only on the respective interpretation of the models. As 
in the case of apparent rotation (§4.1), the tidally induced el- 
lipticity observed by 1PP95I pertains to the specific case of a 
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FIG. 11. — Upper panel: Isodensity contours for a model satellite right FIG. 12. — Number density, velocity dispersion and mean RV trends for the 

before it becomes completely unbound (5 % of initial mass still bound). As satellite shown in the upper panel of Figure 1 1 . Data points are overplotted for 

the satellite is on the verge of being destroyed the bound component of the comparison as open symbols. The model satellite has become significantly 

satellite enlarges and elongates. Lower panel: Isodensity contours for one enlarged compared to its initial size, and it has been scaled down by a factor 

of the Carina's best fitting models (5021). Only particles within two King of four to compare to the data. While the velocity dispersion along the major 

radii, marked by the dotted ellipse) have been considered to mimick the sur- axis of the model seems to provide a good match to the data (after being 

face brightness limitations of current photometric surveys. The contours have scaled down), the rest of the observed properties of Carina are not well match 

been plotted roughly evenly spaced in radius. The density countours show a by the model (even when scaled down), 
remarkably regular structure, with no tail-like structures visible despite the 

satellite losing nearly half its initial mass. That none of our surviving model satellites present tidal 

elongation within their tidal limits suggests then that perhaps 

satellite being observed just before it becomes completely de- the ellipticity observed in some Galactic dSphs corresponds 

stroyed and when the MW has been modeled as a point mass. to a process other than tidal stripping, even if the satellites 

As presented in §3.2.4, model satellites in our survey that are are currently experiencing mass loss. For example, if dSph 

near complete disruption, with only a small fraction of the morphologies ar e the result of tidal s ti rring o f disky dwarf ir- 

initial mass still bound (typically 5% or less), also "puff up" regular galaxies (Ma yer et al]l2001albl. 120071) the triaxiality of 

relative to their original size and they also typically elongate the progenitor system may be preserved for some time. This 

along their major axis (Figure 11, upper panel). might also explain why, in the cases where proper motions 
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FIG. 13. — Isodensity contours for the same model shown in the lower 
panel of Figure 1 1 , but when particles over a four degree radius are consid- 
ered. Again, the contours have been plotted roughly evenly spaced in radius 
along the minor axis. The contour levels shown (from faintest to brightest 
up to r;, m ) correspond to the following surface brightness levels: 35.5, 34.9, 
34.6, 34.2, 33.8, 33.3, 32.0 and 30.4 mag arcsec - ' (assuming the central sur- 
face brightness of the model satellite is Eo=25.5 mag arcsec - ' in the V-band). 
This time, the development of tidal tails is clearly observed, but over scales 
much larger than the King limiting radius. To illustrate how far out one needs 
to go to detect clear tidal tails, current Carina RV members have been over- 
plotted (as solid circles) for comparison. 

have been measured, there is not always a correlation between 
the observed posit ion angle of the satellites and t he direction 
of their orbit (e.g., iPiatek et al.ll200l 120051 120061) . 

Once again, our more systematic exploration of parameter 
space illustrates that one must be careful not to extrapolate 
the results of the earlier, more specific models and expect ef- 
fects that are not necessarily general in nature or that are actu- 
ally beyond the range or sensitivity of current surveys. In this 
case, the elongation reported by PP95 and in our simulations 
is from the nascent tidal tails, and therefore, to see the effect, 
one must look at radii where the debris makes a substantial 
relative contribution to the bound particles (e.g., Fig. 13). This 
is beyond the typical surface brightness limits and/or radii of 
most previous photometric studies of d Sphs (other than our 
own — e.g. , Majewski et al. 2003, Mu noz et al.[ 2006a. 



iSohn et alJl2007h 

4.3. Misperception 3: No Detected Tidal Tails Means No 
Tidal Disruption 

One of the most unmistakable signatures of a system un- 
dergoing tidal stripping is the development of debris tails 
which, in gener al, closely trace the orbit of the satellite (e.g., 
lJohnstonlll998l) . Perhaps the most striking examples in the 
MW system of satellites exhibitin g tidal tails are the Sagit- 
tarius dSph (Majewski et al. 2003) an d the Palomar 5 glob- 
ular cl uster (lOdenkirchen et all 120031: iGrillmair & Dionatosl 
l2006al) . Several other recently found stellar streams have no 
clear progenit or but have been traced for tens of degrees on 
the sky (e.g., IG rillmair & Dionatosl l2006bl: iGrillmairi 120061: 
iBelokurov et al.l 120071) . Several photometric studies have 
been carried out over the years to look for unmistakable ev- 



idence of tail-like features around other Galactic dSphs, bu t 
with the possible exception of Carina dMunoz et al.l l2006a). 
such evidence has been lacking. This absence of detectable 
tail-like structural perturbations at large radii, especially when 
combined with a lack of the (inappropriately) expected ro- 
tation signature of PP95, has often prompted the conclu- 
sion that Galactic tides are unimportant for most dSphs (e.g., 
Walcher et alj|2003l IColeman, Da Costa. & Bland-Hawthornl 
20051: iSegall et al.l 120071) . However, we can use the results 
of our extensive A^-body simulation survey to investigate the 
feasibility of actually detecting tidal tails around satellites un- 
dergoing tidal disruption. 

To clarify the expectations of finding tidal tails we simulate 
"observations" of our A^-body models (where we know with 
certainty that tidal tails are present). If one assumes that the 
central surface brightness of our best matching models (5021 
and 502G) corresponds to the observed cu rrent value f or Ca- 
rina (E =25.5 mag arcsec" 1 in the V-band. iMateoll 1998b . then 
the surface brightness of their density at r/ lm would be £ > 
31 mag arcmin" 1 , which is beyond the sensitivity of typical 
photometric methods used to explore the surfac e brightness 
distribution of dSphs (see iMajewski et aT1l2005l for detailed 
comparisons of these methods). Indeed, in their exploration 
of the forma tion of a MW-like galaxy h alo by A-CDM infall 
of satellites, iBullock & Johnst on (2005) found most of their 
resulting debris streams to be fainter than 30 mag arcsec" 2 , 
and K. V. Johnston et al. (in preparation) predict there to 
be only one stream brighter than this on average . The fact 
that our own surveys (e.g., Maje wski et al.l l2000l et. seq.) 
can probe these and the even fainter surface brightnesses of 
tidal tails is because of the highly tuned selection for dSph 
giant stars, which yields much higher signal-to-background 
sensitivity in the very diffuse parts of the dSph systems. Fig- 
ure 11 (lower panel) shows that "observations" of one of our 
best matching Carina simulations (model 5021)), to a surface 
brightness comparable to those of most current dSph photo- 
metric surveys would not have detected the presence of tidal 
tails around the model satellite. 

Figure 13 expands by a factor of four the spatial scale of the 
density contours shown in the bottom panel of Figure 1 1 and 
includes even lower surface brightness features to the point 
where tidal tails are now visible (the dotted line delineates 
r/ lm for this model, a radius that is well within the regions 
where there are tails). This figure demonstrates another chal- 
lenge that photometric searches for tidal tails face: the de- 
bris around the bound core of the dSph does not organize 
into coherent tails immediately outside r/,„, but the presence 
of tidal tails is only evident over physical scales much larger 
than this. Therefore, much larger systematic photometric sur- 
veys to very low detection thresholds are needed to detect un- 
mistakably obvious tidal tails from real dSphs. On the other 
hand, unbound tidal debris does collect more or less radially 
symmetrically surround the dSph just beyond the tidal radius, 
and it is the presence of this "extended population" structure 
in the density profiles of satellites that is the more obvious, 
and more likely to be detected, expected photometric signa- 
ture of tidal disruption. As pointed out in Section 1.1, most of 
the classical MW dSph satellites do exhibit this very feature 
in their density profiles. 

4.4. Misperception 4: No S-shape Morphology Means no 
Tidal Disruption 

Another common expectation for tidally disrupting sys- 
tems related to the presence of tidal tails is the de- 
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FIG. 14. — Isodensity contours for the same model and particles shown in 
Figure 13, but this time from a line of sight nearly perpendicular to the plane 
of the orbit. As it can be seen, the expected, tidally induced .S-shape is barely 
discernible. 

tection of stellar density contours that become progres- 
sively more elliptical with radius along with a varia- 
tion of their position angles, which together result in 
the typical S-shape morphology associated with tidall 
stripped systems (e.g.,|PP 95: Comb es. Leon. & Meyla nTl99 
IChoi. Guhathakurta. & Johnston! I2002T) . Examples of these 
morphologies have been detected in Local Group systems al- 
ready proy^njobeisnrptingj^^ 1 satellite NGC 
205 (IChoi. Guhathakurta. & Jo hnston 2 0021) and the Galactic 
globular cluster Palomar 5 (Odenki rchen et all2003h . 12 Thus, 
given the apparent connection between S-shaped morpholo- 
gies and tidal disruption, one might be tempted to use the lack 
of detection of this morphology to argue that a system has not 
undergone tidal interaction with our Galaxy. However, closer 
examination of cases where this S-shape morphology has been 
detected shows that usually, with the clear exception of Ursa 
Minor (discussed further below), the isophotal twisting occurs 
over physical scales larger than r;,,„. 

In addition, previo us simulations (e.g. 

ICombes. Leon. & Meyla nfl999l) as well as our own presented 
here demonstrate that, depending on the orientation of the 
satellite and its tails to our line of sight, this S-shape mor- 
phology may not be seen at all. Our best matching models 
for Carina demostrate this situation. In Figure 13, where the 
tidal tails of the satellite are clearly seen, no classical S-shape 
is observed. In this case, the lack of detection is readily 
explained by the fact that this isophotal twisting occurs in the 
plane of the orbit, and, given the Galactocentric distance of 
Carina, which like most of the classical Galactic dSphs, is 
large (> 60 kpc for most cases), the projection of the orbit 
on the sky is not favorable for observing this morphological 
feature (i.e., the Sun lies relatively close to the orbital plane 



for distant Galactic satellites). Moreover, even if viewed in 
the correct orientation, the S-shape might be very subtle. 
Figure 14 shows the density contours for the same model and 
particles shown in Figure 13, but this time from a line of sight 
nearly perpendicular to the plane of the orbit (e.g., projected 
onto the Galactic Yqq-^cc plane). Even in this preferred 
orientation for observing the S-shape it is barely visible. 

T he case of Ursa Mi nor, which shows an S-shape morphol- 
ogy dPalma et alJ2003h but where the reported S-shape is seen 
well inside its r/; m , is thus unique and i ntriguing. Th is dSph 
is at a heliocentric distance of 66 kpc (Mateo| [l998l) . and so 
it is unlikely that its orbital projection on the sky is favorable 
for viewing the S-shape, regardless of the shape of the orbit. 
Moreover, even if this were the case, we would expect to see 
isophotal twisting over scales larger than r/,- m . In this particu- 
lar case, the S-shape may be a coincidence, with the isophote 
bending the result of other physical effects at work to produce 
internal substructure. For example, the well-known secondary 
peak in the 2D stellar density distribution of Ursa Minor, dis- 
placed ~ 300 pc from its cen ter, and shown to be dynamically 
cold by iKlevna et al.l d2003f) . could potentially result in iso- 
pleths that just happen to mimic the S-shape twisting even if 
they are not of tidal origin, at least in the classical way. 

4.5. Misperception 5: Flat/Rising Velocity Dispersion 
Profiles Require Extended Dark Halos 

Most current RV surveys of dSph satellites yield veloc- 
ity dispersion profiles that are more or less flat or ris- 
ing all the way to the r/„„ of the system. Recently, 
these velocity dispersion trends have been attributed to 
non-MFL, multi-component mass distributions in dSphs, 
with perhaps some possi ble additiona l contribution from 
velocity anisotropy (e.g., [Eo kas 2002; Klevna et al.j 120021 : 
[ Lokas. Mamo"n& Prada 2005; Mashche nko et al.l2005Ll2006t 
I Walker et al.l l2006at iKoch et al.l l2007allbl) . These multi- 
component models yield both total dSph masses that are up 
to an order of magnitude higher than those obtained by tra- 
ditional core fitting techniques, as well as physical scales 
much larger than those of the observed light distribution, 
suggesting that dSphs are currently surrounded by extended 
dark matter halos. In some cases the derive d satellite to- 
tal masses approach the 10 9 - 10 10 M M range dLokasI 120021; 



12 Other s ystems, like the An d I (McConnachie & Irwin 2006) and Ursa 
Minor dSphs I Palma et al. 2003) show this S-shape morphology, but in these 
cases, the presence of this feature has been adopted as one of the arguments 
that these dsphs are tidally disrupting. 



lMashchenkoeiini200ll2006l:lRead et alJl2006b 

However, our analysis of dSph properties within the context 
of MFL, tidally disrupting admits a workable alternative inter- 
pretation of the velocity dispersion profiles of dSph galaxies. 
We summarize our results with respect to the behavior of the 
velocity dispersion profile as follows: 

• The most important observable dynamical effect of 
Galactic tides is to inflate the velocity dispersion at 
large radii and to produce a velocity dispersion profile 
that can go from slowly declining to flat or even rising 
depending on the amount of tidal debris released and 
disruption history of the satellite. 

• The central velocity dispersion of systems that retain a 
bound core remains unaffected by tides, which means 
that its value is a good reflection of the underlying 
bound mass of the satellite, even when substantial tidal 
stripping may be ongoing. 

• Completely flat velocity dispersion profiles are, in gen- 
eral, only achieved when the satellite has become com- 
pletely unbound. The size of the velocity dispersion is 
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driven by a complex combination of the initial satellite 
mass, density, orbit and number of prior orbital pas- 
sages. 

• The inflation of the velocity dispersion as a func- 
tion of radial distance is not necessarily an effect 
of "cylindrically-averaged" (i.e., azimu t hally av eraged) 
dispersions as suggested bv lRead et alj d2006al) . Intrin- 
sic inflation is observed even when the velocity disper- 
sion is studied along the major axis (see Figures 3 and 
6). 

• Completely or mostly bound systems show velocity dis- 
persion profiles that peak at the center and decrease 
with radius. 

Thus, there are at least two viable explanations for the 
flat/rising velocity dispersion profiles of dSph galaxies. How- 
ever, beyond just explaining the velocity dispersion profiles, 
our MFL tidal disruption models are the first models that try 
comprehensively to match all presently observed general fea- 
tures of dSphs (except ellipticity; see §3.2.4), and that satisfy 
all specific constraints of the presently best-studied dSph in 
particular. No multi-component, tidal-stripping-free model- 
ing of the Carina dSph has yet been attempted (several studies 
have tried to match other dSphs with more limited data sets, 
but Carina is currently the system with the most extended ra- 
dial coverage, besides the Sgr dSph), so it remains to be seen 
whether such models can successfully match these observa- 
tions as closely and with as much simple elegance as tidal 
disruption models. 13 

4.6. Misperception 6: Dark Matter and Tidal Disruption in 
dSphs are Mutually Exclusive 

Previous investigations of tidal effects in dSphs have of- 
ten proceeded with an "all or nothing" approach with regards 
to the dark matter content of these systems. For example, 
iKuhnl (119931). iFleck & Kuhnl (120031). fPP95l IOLA95L iK rouna 
d 1997b and lGomez-Flechoso & Martmez-Delgadol d2003l) . all 
explored dark-matter-less systems, while most of the refer- 
ences in the previous subsection generally sought to explain 
the outer properties of dSphs (most often the dynamical, but 
not necessarily the density, properties) entirely with extended 
dark matter halos. Our approach has been to explore dark- 
matter-dominated systems, but ones where tidal disruption is 
still a prominent feature shaping their dynamical and struc- 
tural character. 

Regarding the mass content and M/L of our test case, the 
Carina dSph, our modeling yields revealing insights: 

• Carinas with no dark matter (current M/L ~ 1 — 3, or 
M ~ 1 x 10 6 M ) are hard to make: Low mass mod- 
els on benign (rather circular) orbits with scale lengths 
comparable to real dSphs and that manage to survive 
(for 10 Gyr or five orbits) do not show the proper den- 
sity nor velocity dispersion profile. On the other hand, 
low mass models on highly radial orbits are gener- 
ally destroyed relatively easily (as has previously been 
found by PP95) and in order to make such satellites sur- 
vive until the present time the initial densities have to be 
so high that the scale lengths are unrealistically small 

13 IWJt2b07l) has created multicomponent dSph models that produce both 
flat velocity dispersion profiles and breaks in the density profiles of several 
other dSphs with less radially extensive velocity data. 



compared to dSphs. Moreover, during those evolution- 
ary phases when the models still retain a bound core, 
the central velocity dispersion is also too low to provide 
an acceptable match to the data even if one resorts to the 
effects of binary stars or atmospheric jitter in the giant 
star atmospheres to account for the differences between 
the model and observed central velocity dispersion (see 
§1.1). 

• On the other hand, MFL Carinas with current masses 
significantly higher than 3 - 4x 10 7 M Q are impossible 
to make conform to the observed velocity dispersion 
and density profiles simultaneously. They are harder to 
disrupt compared to lower mass satellites for a given 
orbit and length scale, and even when they do disrupt, 
the debris is often too energetic compared to the ob- 
servations (§2.2). In these cases, if one tries to repro- 
duce the Carina number density profile (by varying the 
mass loss rate) the models need to be put in very eccen- 
tric orbits, which, in turn, "overheat" the particles and 
produces velocity dispersion profiles that significantly 
depart from that actually seen in Carina. Equivalently, 
if one tries to match the velocity dispersion profile by 
increasing the density (and therefore produce a modest 
mass loss rate), then the physical scale of the satellites 
becomes unrealistically small. 

• Our best matching models (those with initial masses 
of 3.6-3.8 x 10 7 Mq) independently predict present 
masses for Carina consistent with those derived for 
the r eal dSph from the core-fitting method ( ~ 2 x 10 7 
M g . iMateo etaDll993h iMufioz et al.ll2006al) . At first 
this may seem unsurprising, given that the initial mass 
and length scale of the model Carina were specifically 
tuned so that the final stage of the simulated satellite 
matched its observed central velocity dispersion and 
core radius. It might be argued that the agreement 
between the estimated mass of Carina (based on its 
central velocity dispersion) and that of the simulated 
model is just a reflection of the fact that our models 
meet (by design) the assumptions used by the core- 
fitting technique: a single-component, isotropic, MFL, 
spherically-symmetric distribution in equilib rium fol- 
lows the distribution assumed bv llllingworthl i 19761) for 
estimating the total mass. 

But this agreement between model and observations is 
not trivial. While providing by design a good match 
to the observed central velocity dispersion and scale 
length, the masses for the best-fitting model Carinas 
also provide an excellent match to the velocity and 
velocity dispersion profiles at large radii (when rep- 
resented either as a function of radial distance or po- 
sition along the major axis) out to the limits of the 
current observations. In addition, the amount of de- 
bris released, which is also a function of the instan- 
taneous mass content, provides a very good match to 
the observed number density distribution. Interestingly, 
the core fitting technique, applied to the model Carina, 
yields the proper instantaneous bound mass value even 
after the satellite has lost most of its initial mass, and 
this technique for gauging the mass only breaks down 
at the point of near total destruction of the satellite. At 
all intermediate time steps (as Tables 3 and 4 show) the 
central velocity dispersion effectively tracks the bound 
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mass content of the satellite. 

In the end, the ability of our best fitting model to simultane- 
ously match all the observed properties of Carina (except, of 
course, th e observe d ellipticity, §3.2.4) leads us to the same 
finding as IOLA95I in that the central velocity dispersion of 
modeled dSphs generally reflects their equilibrium value even 
if the satellite is undergoing significant tidal disruption; there- 
fore core-fitting is validated as a useful mass indicator if, even 
in the presence of tidal disruption, dSphs are adequately rep- 
resented by MFL models. In the case of Carina, we show that 
the dSph is well described by tidally disrupting MFL models, 
even though it is still a very dark matter-dominated satellite. 
Clearly, it is possible for dark matter-dominated models to ex- 
hibit luminous tidal disruption. 

4.7. Misperception 7: The King Limiting Radius is 
Meaningless as a Tidal Radius Indicator 

As discussed in §4.7, at present, the underlying dark matter 
distribution in dSphs is still unknown, and we are left only 
with the distribution of the luminous matter to understand the 
structural properties of these systems. A relevant question 
for the MFL models explored in this paper then is, if in fact 
present dSphs can be modeled well as single-component ob- 
jects, does rn m provide any useful information regarding the 
true tidal radius of the satellite? 

The tidal radius of a satellite can be defined as the place at 
which stars within the satellite become unbound to the system 
and, in turn, become bound to the host galaxy. Ear l y studies 
of glo bular clusters dvon Hoerneii[T957l:)Kinglil962l:IOh et alj 
1 19921) derived simple relations to estimate this boundary that 
depended only on the mass of both the satellite and the host 
galaxy as well as the shape of the satellite orbit, but soon 
afterward it was appreciated that the tidal radius of a sys- 
tem was also a function of the type of orbit the stars within 
the satellite have; stars on prograde orbits are more easily 
stripped than stars on radial orbits, and the latter are more eas- 
ily stripped than stars in retrogr ade orbits dToomre & Toomrel 
1972t [Keenan & Innanen|[T975l and more recently iRead et all 
2006a). Thus, it became clear that the true tidal radii of 
dSphs are not a well defined edge, but rather a shell-like zone . 
Previous /V-body simula tion studies (e.g., iJohnstonl Il998t 
iKlimentowski et aTl2 007) and the simulations presented here, 
show that unbound debris can linger deep inside the space de- 
fined by the bound population and, at the same time, stars still 
bound to the satellite can remain well into the region domi- 
nated by tidal tails, effectively "blurring" the "tidal boundary" 
from a Roche-like surface to a population overlap zone. 

Despite the impossibility to define a clear, well delineated 
tidal boundary, if one knows the bound mass for a satellite in 
a given orbit with eccentricity e, the tidal rad ius at perigalac- 
ticon can be estimated using dOh et alJ [f992). 
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where a is the orbital semimajor axis for the satellite, and 
MdSph and Mq are the mass of the dSph and the MW (inside 
a), respectively. The question is how does this compare to the 
K ing limit i ng rad ius ? 

Burkert (1997) tried to estimate the enclosed mass, Mds P h, 
of several dSphs (including Carina) using this equation and 
adopting r;„„ as R t idai- He then compared the derived Mjsph to 



the dynamical masses (those derived from core-fitting), and 
argued (1) that the latter masses are much too high compared 
to the former masses, or (2) reversing the problem, if the 
masses of dSphs are indeed similar to their dynamical masses, 
then R/idai > > ru m and therefore r;, m has no meaning as a tidal 
boundary. In the latter scenario, no extratidal stars should be 
observed in dSphs because the luminous component is effec- 
ti vely shie lded insid e a much larger tidal boundary. 14 

IOLA95I and IPP95I also explored the relation between ru m 
and R, 

idcd an d found that a satellite is able to survive tidal in- 
teraction (i.e., not be completely destroyed by tides) for sev- 
eral Gyrs only if its r/,- m is smaller than twice its R t idai- 

The results from our survey provide us with a possible way 
to weigh in on the question of R,id a i versus r?, m . In Tables 3 
and 4, we show the instantaneous tidal radii of our satellites 
as determined by equation (4) using the bound mass of the 
system at that point of evolution and the density of the MW 
at thos e apogalactica. As may be seen, our models meet the 
IOLA95I conditions for luminous tidal stripping (r;, m > Rudai) 
and survival (2R,jd a i > ru m ) at the same tim e. However, ou r 
results are at odds with the conclusions of Burkert (1997): 
Despite having masses of a few times 10 7 (as correctly in- 
ferred by core-fitting in our models), our best fitting models 
have a tidal radius smaller than its rn m and do show extended, 
unbound populations. This discrepant result is explained by 
an arithmetica l mistake in the application of equation (4) by 
lBurkertl(fT997l) . 15 

If single-component models undergoing tidal stripping are 
adequate enough that they provide a functional description 
to the observed properties of dSphs, then they revalidate the 
value of r/,„, as at least a crude approximation for an upper 
limit for the R ida i of dSph systems (i .e. 2R, ida i > r Um >R ti dai)- 
In this regard, Maje wski et al.l (120031) have already argued that 
in the case of Sagittarius, the only dSph clearly in the process 
of being tidally disrupted, the true tidal radius has to be well 
inside r/„„ of the light distribution. 

4.8. Misperception 8: The Combined Observations ofdSphs 
Show That They Are Not Being Tidally Disrupted 

The main purpose of this work has been to demostrate that 
Galactic tidal effects can produce the observed properties of 
dSphs, in particular, those of Carina, the Galactic dSph with 
the most extensive structural and kinematical data sampling 
to date, but otherwise a system not unusual compared to most 
other dSphs in its structure and dynamics. We have also pre- 
viously reported tha t the Leo I dSp h is well-fitted by tidal 
disruption models dSohn et al.ll2007l) . and possesses charac- 
teristics (an asymmetric velocity distribution at large radius) 
that cann ot be readily explained by an exten ded dark halo. 
Recently, Mat eo. Olszewski. & Walkerl (120071) reported a ro- 
tational signal in the outer parts of Leo I they interpret as the 
signature of tidal stripping, in line with our arguments for Leo 
I. While the Sgr dSph is an uncontroversial case of a dSph in 
the p rocess of being accreted, ther e is a tendency to attribute 
Sgr (lMateolll998l iGilmord 12004). and now Carina as well 
(Gilm oreet alJl2007lT as "exceptions to the rule" that dSphs 

14 This does not preclude dSphs from having extended luminous compo- 
nents if these form part of an extended, bound halo. 

15 Whi le rewrit i ng e quation (4) to solve forMo, which involves cubing the 
equation, Burkert 1 1997) did not cube the expression within the braces on the 
right-hand side where the effect of the orbital eccentricity e is included. The 
error becomes increasingly significant as the orbital eccentricity of a system 
increases, but even the value of the expression for the case where e = is 
incorrectly reported in that paper. 
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FIG. 15. — (a) Number density profile, (b) velocity dispersion versus angu- 
lar distance and (c) mean RV trend along the major axis for model 5021 when 
only particles within rn m are considered. This mimics the typical surface 
brightness limits of "traditional" dSph photometric + spectroscopic studies. 
The dotted line i n panel (a) repr esents a typical background level from pho- 
tometric surveys iMaiewski et al. 200^). 



are in general not disrupting. The results presented here sug- 
gest that it is not straightforward to argue against this position 
with current data on dSphs other than Sgr because, even if 
dSphs are undergoing significant tidal stripping, the effects 
can be for the most part undetectable (e.g., §4.3). 

To illustrate how the lack of sufficiently sensitive data can 
devolve into an apparently definitive argument against tidal 
disruption we use recent discussions regarding the state of 
the Draco dSph galaxy. The lack of concordance between 
observational results and expectations from previous TV-body 
simulations have been used to conclude that Draco cannot be 
significantly tidally influenced by the MW. Th e main argu- 
ments supporting this view are summarized by iKlevna et alj 
d2001l) . who argue that Draco cannot be significantly affected 
by Galactic tides because: (1) no significant rotation is ob- 
served within the radial extent of the current kinematic al sur- 
veys (which are limited to Draco's r/ ,•„,), (2) Draco shows a 
"roundish" appearance (low eccentricity of ~ 0.3) in appar- 
ent contradiction with th e predict ions of the TV-body simula- 
tions done by IPP95I and IOLA95L and (3) its surface bright- 
ness profile is adequately fitted by a Plummer law without the 



need to resor t to an "extended" populati on in the outer parts. 16 
Furthermore, Oden kirchen et alJ (1200 ll) add that n o tidal tails 
from D raco are visible within the SDSS data and lRead et alj 
(2006b) argue that tidal stripping of dSphs should cause the 
velocity dispersion profile to rise rather than to remain flat, 
which is the observed behavior for Draco. 

While not claiming by any means that Draco is being tidally 
disrupted, we have shown in this paper that an object with 
Draco's observed structure and kinematics can still be experi- 
encing tidal stripping of its luminous matter. To illustrate fur- 
ther how this fact may be disguised, Figure 15 shows the den- 
sity profile, velocity dispersion profile and RV trend (in GSR) 
of one of our best matching model Carinas, but with data "ob- 
served" only out to r/,,„, as is the case for Draco (compare to 
Figure 10, where the full model data are shown). 17 Carina and 
Draco are at nearly similar distances (80 vs 100 kpc), have 
similar sizes, similar luminosities and both have high central 
M/L. However the two dSphs have different vgsr, with the 
velocity suggesting a vastly larger apogalacticon for Draco; 
nevertheless our Carina model is useful for this illustrative 
exercise. 

In Figure 15a we show that, as observed in the case of 
Draco, a Plummer model does an excellent job fitting the pho- 
tometry, even though, in fact, the outer parts of the model den- 
sity profile are dominated by debris and the bound component 
has departed from the initial Plummer distribution. Figure 
15b shows that the velocity dispersion profile does not signif- 
icantly increase beyond the central value within the radial ex- 
tent of our "observ ations", a result that contradicts the above 
iRead et all (l2006bl) claim for the expected rising behavior of 
the dispersion profile used to argue against tidal disruption in 
Draco. Figure 15c shows that no statistically significant ro- 
tation is observed within the nominally probed Draco radius, 
where signi ficant means greater than the velocity dispersion 
dPrvorlll996l) . Finally, the intrinsic eccentricity of the bound 
component of the model, as "observed" in Figure 11, is close 
to zero. Both of these latter findings are also contrary to ar- 
guments by Kleyna et al. (2001) against Draco disruption, 
ar guments bas ed on e xpectations from the models developed 
bv lOLA95l and lPP9l 

Taken together, the Figure 15 "observations" of a tidally 
disrupting (fractional mass loss rate of 0.1 per Gyr), MFL 
satellite model — when coupled with common mispercep- 
tions regarding the expected characteristics of such struc- 
tures — can be used to argue compellingly that our model 
satellites are neither tidally disrupting nor well-expained by 
MFL models. With this apparent " failure of the tidal model" 
(E 

essen. Grebel. & Harbeck 2003J), it is understandable how 
searches for an alternative mechanism to explain the flat 
Draco velocity dispersion profile would be deemed unneces- 
sary or irrelevant. 

5. FINAL REMARKS 

5.1. What Other Studies Are Saying 

As mentioned in the Introduction, the IOLA95I and IPP95 
studies provided important benchmarks in the exploration of 
the structure and dynamics of dSph galaxies, but in many 



16 More recently, Wilkinson et al. ( 2004) found that the surface brightness 
profile of Draco seem s to show a "break " at about 25', which is inside its 
observed r/,„,, whereas Segall et al. 1 2007) could not confirm the presence of 
this extended population down to three orders of magnitude below the central 
number density. 

17 The kinematical data on Draco in fact barely reaches its r;,-,„. 



22 



Modeling Dwarf Spheroidals with Dark Matter and Tides 



ways the results of these studies have been taken out of con- 
text and extrapolated to regimes that, we believe, were not 
actually intended by these authors. The primary emphasis of 
these two studies was to establish whether tides could inflate 
the central velocity dispersions of dSph galaxies as a way to 
eliminate the necessity of dark matter; the goal of these stud- 
ies was not to rule out that any tidal effects could be occurring 
in dSph satellites. Even though PP95j and lOLA95l succeeded 
in demonstrating that in fact tidal effects don 't typically in- 
flate central velocity dispersions (a result we have also found 
here), it doesn't follow that tidal disruption is not happening 
in dSph systems. Nor is it appropriate to take the observ ed ef- 
fects in the v ery well-defined and restricted range ofthe |PP95l 
and IOLA95I simulations and apply them generally to other 
contexts. 

For example, perhaps the most cited broadside used to sink 
the tidal disruption pr oposition is the lac k of observed dSph 
rotation, — e.g., from lKoch et al.1 (l2007bh : "... the most effi- 
cient way of gathering evidence whether a particular dSph is 
being affected by tides sufficiently strongly that its kinematics 
are being modified is to look for any sign of apparent rotation 
in its outer regions." This leads to arguments against tides not 
only in particular dSphs — e.g., "all in all, there is no kine- 
matical evidence of any significant velocity gradients in Leo 
II, either due to rotational support or produced by Galactic 
tides. It therefore appears that tides most likely have not af- 
fected Leo II to any significant degree... " dKoch et al.ll2007bh 

— but against tides generally in Local Group dwarf galaxies, 
since "...to date, all but one of the dSphs of the Local Group 
show no significant velocity gradients." But our analysis of 
this manifestation of tidal disruption (§3.1.2 and §4.1) shows 
clearly that rotation is not a reliable test of tidal disruption — 
a r esult th at, in hindsight, should not be that surprising since 
(1) IPP951 already had demonstrated that the effect was pro- 
nounced in only a limited set of circumstances (e.g., at peri- 
galacticon, or when the MW is modeled as a point mass), (2) 
both lPP95l and lOLA95l stress that this phenomenon would be 
visible only at very large radii (i.e., outside the tidal radius, 
typically on > 1 deg scales) where few dSphs have been ex- 
plored, and (3) with small amplitude in most cases. 

Meanwhile, in a ddition to the rotation argument, 
iGilmore et alJ (120071) suggests that flat dispersion pro- 
files also argue against tidal disruption: "...there is some 
evidence ... for flat velocity dispersion pro files to large 
distances (ISohn et alJl2007t iMufioz et al.ll2005h . ... which are 
inconsistent with simple tidal disruption effects, particularly 
since most dSphs show no evidence of apparent rotation 
dKoch et al.l |2007a)." But of course, in addition to our 
modeling here, tidal disruption models have already been 
demonstrated to be capable of producing flat dis persion 
profiles in a n umber of studies in preceding years (PP95; 
iKroupal fT997t iMayer et all 120021) — including in the very 
Sohn et al. study referenced in the above passage. Indeed, 
the dominating influence of unbound stars on the velocity 
dispersion at large radius is one of the very reasons MFL tidal 
disrupti on models have been invoked (a t least by our group 

— e.g. JLaw et al1l2^0llSolmiraT]l2007h as a potential dSph 
scenario. In actuality, it is only in this circumst ance that such 
model s could work, for, as summarized by IGilmore et al.1 
d2007l) : "More generally, there are unavoidable consistency 
requirements in any mass follows light model. In any model 
where mass follows light the projected velocity dispersion 
must be maximum at the centre, and then fall monotonically. 



For a well-mixed (star cluster) system, the velocity dispersion 
will decrease by roughly a factor of two over three core radii. 
This is an unavoidable requirement for any mass follows light 
system, and is observed in star clusters ... Such a velocity 
dispersion profile is not required by data in any well studied 
galaxy, however small, further emphasising the intrinsic 
difference between (virial, King) star clusters and galaxies." 
The latter description as written is obviously correct in the 
absence of tidal effects, but it should not be used as an 
argument against MFL models when they are being used to 
explore the very regime — tidal disruption — to which the 
statement is not applicable. 

Perha ps another unfortunate aftermath to the IPP95I and 
IOLA95I papers is that, by their example of exploring com- 
pletely dark-matter-free satellites, these two studies may have 
had an unintended and undue influence in paralyzing sub- 
sequent discussions into considering primarily two polarized 
viewpoints or modeling strategies — namely, those where (1) 
the model satellites have no dark matter (thus requiring some 
reassessment or proposed revision of the physics involved 
with conventional mass determinations to reconcile the DM- 
less hypotheses with the observations), or where (2) dark mat- 
ter is invoked to be responsible for all observed effects seen in 
dSph galaxies (most starkly demonstrated by the frequent in- 
corporation of large dark matter halos to explain the structure 
and dynamics of dSphs at large radii). In our opinion, and as 
demonstrated by the results presented in this paper, there is an 
enormous middle ground between these polar positions. 

Nevertheless, more recently the weight of discussions has 
strongly tacked to the side of the large dark matter halo hy- 
pothesis to explain the structure of dSphs — to the point 
that there are often strongly expressed sentiments that there 
simply is no longer any reason for consideration of alter- 
native models. Ironically, this viewpoint is partly driven 
by, on the one hand, the apparent success of A-CDM in 
explaining structure on large scales, but on the other hand 
the lack of success of A-CDM in explaining structure on 
galaxy scales so that more massive satellites are needed to 
help resolve the missing satellites problem. Extended, non- 
MFL dark matter halos can not only account for some of 
the observed properties of dSph galaxies at large radii and 
fit into A-CDM paradigm, but they mitigate to some de- 
gree the mismatch in the observed and A-CDM-predicted 
mass spectrum of satellite systems. Thus, for example IWul 
d2007l) concludes "Because our [multi-component, extended 
dark matter] models can fit both radial velocity profiles and 
surface number density profiles, the so-called 'extra-tidal ex- 
tensions' in the surface num ber density profiles found by 
llrwin & Hatzidimitrioul d 1 9951) and lWilkinson et all d2004l) do 
not require any special ad hoc explanation. Thus it is not valid 
to consider them as the evid ence of tidal stripping, as pro- 
posed by some authors (e.g., Martmez-Delgado et al.l 120011: 
Gomez-Flechoso & Martmez-Delgadg 120031 : IMufioz et alJ 
20051) .""* 

But we have shown here that our tidally disrupting MFL 
models can fit not only the radial velocity profile and sur- 
face number density profile of the Carina dSph but also the 
velocity dispersion profile. In addition, tidal disruption is a 
naturally occurring, established phenomenon in the evolution 
of galaxies (see §5.2). Thus, the MFL models are less ad 
hoc than implied above and, indeed, they provide a fully self- 
consistent explanation for all of the structure and dynamics 

18 Emphasis added. 
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of at least the Carina dSph. Unfortunately, what such mod- 
els are not consistent with is the prevailing understanding of 
dark matter, wherein the latter is non-baryonic and does not 
strongly couple with luminous matter. Thus it is not clear why 
mass should follow light. 

On the other hand, while multi-component models are 
aligned with the A-CDM orthodoxy, there are, as yet, no spe- 
cific predictions in the context of dSphs from A-CDM the- 
ory for the expected structure and dynamics of the luminous 
matter embedded in dark halos. Indeed, the inability to pre- 
dict both the observed cored profiles and the angular mo- 
mentum distributions of galaxies are ack nowledged additional 
failing s of CDM on galactic scales (e.g jNavarro & Steinmetzl 
l2000allbl) . and multi-component models of dSphs must also 
adopt empirically-based, ad hoc distrib utions for the lumi- 
nous matter, such as King profiles (e.g.. | Strigari et al1l2007h 
iPenarrubia. Navarro. & McConnachidl2007l) . 

Clearly then each the MFL and multi-component models 
have unsatisfying philosophical aspects and neither can ac- 
count for why luminous matter should organize in a particular 
way in the presence of substantial dark matter; under such 
circumstances we believe it is pertinent to address more ba- 
sic questions such as "Which model works better?" or "Can a 
given model be definitively ruled out?". One of the goals of 
this paper is to try to "clear the air" with regard to (MFL) tidal 
disruption models, i.e. to clarify their expected characteristics 
over a wider (and more likely) range of parameter space than 
previously explored, so that, ultimately, stronger tests can be 
developed for either verifying the tidal disruption hypothesis 
or ruling it out. Ironically, while we can make testable predic- 
tions for the tidal disruption hypothesis that can be checked 
with the discovery and observations of dSph stars at ever large 
radii, there is no direct way to check extended dark matter 
models at large dSph radii, beyond where there are sufficient 
luminous tracers and where most of the dark matter halo is 
supposed to lie. In this sense the extended dark matter models 
are virtually unconstrained, because the extent and mass of 
the dark matter halos depends strongly on how one constructs 
the mass function (or v Ci y c distribution function) beyond the 
light. And, even at radii where some stars are present, with 
extra free parameters at their disposal, multi-component mod- 
els can usually be constructed to accommodate the properties 
of those stars. Perhaps this is why tidal disruption models are 
"more vulnerable" in the common discourse on dwarf galax- 
ies: Whereas the challenging exercise of discovering dSph- 
associated stars at large radii is the sine qua non of the tidal 
disruption scenario, extended dark matter halo theories are es- 
sentially immune to whether stars are found at large radius 
or not. At this point it seems that virtually the only way to 
disprove the existence of an enormously extended dark matter 
halo around any particular dSph is to push the discovery of as- 
sociated stars to reductio ad absurdum radial limits where un- 
reasonable dSph masses are needed to kee p the stars bound to 
the sat ellite (a strategy employed in, e.g.. iMunoz et al.ll2005L 
l2006ah . 

5.2. What We Are Saying and What We Are Not Saying 

I n this pap er we attempted to go beyond the seminal PP95 
and IOLA95I papers to explore tidal disruption models across 
a wider, and we believe more realistic, part of the parameter 
space of satellite mass, density, and orbit. Like these land- 
mark papers, we used one component models, but with a key 
difference in that we adopt satellite masses consistent with 
dSphs having significant dark matter contents — i.e. consis- 



tent with the central mass-to-light ratios. 

Remarkably, despite the simplicity of our models, which 
by virtue of having a single dynamical component invoke the 
minimum number of free parameters (satellite mass, linear 
scale, and orbit), we have been able to zero in on examples 
that rather faithfully reproduce Carina's velocity profile, ve- 
locity dispersion profile, density distribution and other prop- 
erties over its entire sampled radius, and without appealing 
to special timing or other particular viewing conditions (ex- 
cept that the satellites match the current properties of the real 
Carina system). We have previously also successfully mod- 
eled the pr operties of the Le o I dSph system using similar 
techniques (So hn et afll2007l) . 19 In both cases, the defining 
characteristic of the successful family of models is high or- 
bital eccentricity. And in both cases the central velocity dis- 
persion is found to faithfully represent the bound mass and 
central M/L of the parent satellite (40 and 5, respectively, 
for Carina and Leo I), while the properties of the outer part 
of the satellite (the density law, velocity profile, and velocity 
dispersion profile) are strongly influenced by the presence of 
unbound tidal debris. But, fa r from being ad hoc, as charac- 
terized and dismissed bv lWul (12007). this tidal debris is a nat- 
ural consequence of the dynamical evolution of the satellite 
and, moroever, it represents a phenomenon already observed 
in Nature: Indeed, were one to subscribe to the Principle 
of Parsimony, a successful one component model would be 
preferred over an equally successful multi-component model, 
which could be argued to contain superfluous elements 20 . 

Our own position regarding the implications of the simula- 
tions presented in this paper takes into account not only the 
practical success of tidally disrupting, MFL models in de- 
scribing at least some specific MW dSph systems, but several 
additional, more general observations: 

1. Both dark matter and tides exist in the Uni- 
verse . Tidal effects are expected in CDM mod- 
els (e.g.. lHavashi et alj 120031: iKazantzidis et alj|2004t 
Kravtsov, Gnedin & Klvpin 2004; Bullock & Johnston 
20051: iKnebe et alJl200(fr 

2. The degree of importance of tides must surely vary 
among Galactic satellites, being a function of the struc- 
ture of the satellites and the nature of their orbits. 

3. We see tidal debris streams throughout the Galac- 
tic halo, some of them massive or dynamically 
warm enough t o have come from t h e disruption of 
dwarf galaxies (jNewberg et alj l2002t iMaiewski et alj 
2004: iDuffau et alj|2006t iGrillmair & Dionatosl l2006b; 
Grillmairll2006HBelokurov et al.ll2007l) . 

4. There are growing suggestions of apparent tidal ef- 
fects among MW dSphs besides Sagittarius, Ca- 
rina an d Leo I, including most rec ently in the Ursa 
Major dWillman et alj 12001 l2006h . Ursa Major II 

19 We have used these techniques to model the Sgr dSph system 4Law et al] 
2005) as well, but in this case the focus was less on matching the properties 
of the parent satellite core and more on adequately explaining the properties 
of the tidal debris. 

20 One of the many variants of the philosophical maxim often called Oc- 
cam's Razor is by Thomas Aquinas: "If a thing can be done adequately by 
means of one, it is superfluous to do it by means of several; for we observe 
that nature does not employ two instruments where one suffices." (Summa 
Contra Gentiles, from the Basic Writings of St. Thomas Aquinas, translated 
by A.C. Pegis, Random House, New York, 1945, p. 129). 
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dZucker et al.ll2006bt IFellhauer e t alj |2007l). Canes Ve- 
natici dZucker et al.l l2006al). Bootes jBelokurov et alJ 
2006), and Hercules (IColeman et al.ll2007b systems. 

A theory for the structure and evolution of dSph galaxies that 
can also account for the above, generally non-controversial 
elements in a natural way would have the appeal of an overall 
economy of hypothesis. 

While we have not ruled out large dark matter models 
and our A^-body simulations are by no means intended to be 
unique descriptions of the Carina dSph in particular, or dSphs 
in general, we hypothesize that at least some present dSphs 
might actually be well represented by tidally disrupting, MFL 
models. This hypothesis is not incompatible with A-CDM, 
and it is not our intention to argue against A-CDM or extended 
dark matter halos in general. It would be foolish to do so, for 
it is virtually impossible to prove that extended dark halos do 
not exist, whereas, as we have discussed in §4, it is very diffi- 
cult to prove that tidal disruption is happening in most dSphs. 
Under these circumstances, the most conservative position is 
to say that both possibilities might coexist. 

For example, if most luminous galaxies start out embedded 
in extended dark subhalos of various relative sizes, eventu- 
ally those on more disruptive orbits and/or with lower den- 
sities and/or radial extents will have their dark matter cloaks 
whittled down by tidal mass loss until their luminous matter 
is "exposed" and vulnerable to the same tidal effects (creat- 
ing vis ible tidal debris streams). Indeed. iKlimentowski et alJ 
(2007) show that cosmologically-driven galaxy formation 
models with two-component satellites that also experience 
tidal disruption naturally lead to the evolution of MFL-like 
satellites after their extended dark matter halos have been 
stripped down to the luminous cores, and at the present time 
a certain fraction of Galactic subhalos can be expected to 
be in this state. Thus, the Carina dSph didn't necessar- 
ily start out as an MFL-like satellite, even if it might be- 
have as one now. Other dSphs, such as Draco, might still 
be firmly enshrouded within exte n ded dark matter cocoons 
(lLokasI 120021: iKlevna etal] l2002t iMashchenko et alJ 120051: 
Segall et al. 2007; Wu 2QQ3), although, again, as pointed out 
in §4.9 and Figure 15, were Draco in a tidally disrupting 
phase it co uld be extremely difficult to prove. dSphs like 
Ursa Minor (iMartfnez-Delgado et al.l200UlPalma et alJ2003t 



iMunoz et ai1l2005l) and Leo I dSohn et alj|2007l) might be in 
states similar to Carina as we have modeled it here, whereas 
other systems, such as Sgr and perhaps some of the newly 
found, extremely low luminosity dSphs with evidence of tidal 
disruption (e.g., Ursa Major I and II, Canes Venatici, Bootes, 
Hercules) may be fa rther along this path of metamorphosis. 
IFellhauer et alJd2007l) have already shown explicitly that Ursa 
Major II — the appare nt progenitor of the "Orphan" stream 
(Belokurov et al.l 120071) — is better described by a single- 
component than a two-component structural model, in keep- 
ing with the general scenario outlined here. Finally, we might 
expect there to be some systems now completely destroyed, 
and visible only by their tidal debris. Such debris might be 
in well-correlated tidal streams, or more diffuse "clouds", de- 
pending on the details of the satellite-MW interaction. Indeed, 
we have produced some A^-body models for systems that are 
of low mass, low density satellite on highly radial orbits that 
are completely destroyed, and the remains of these systems 
resemble the large, diffuse, halo "st ar cloud" recently discov- 
ere d in Triangulum- Andro meda by iRocha-Pinto etall (2004) 
and lMaiewski et akl ([2004). Other features, like the "shelves" 
of stars see around the disk of M31, can be cre ated as tidal 
shells by very intense sa tellite-host interactions (Fard al et alJ 
120061: iGilbert etal]l2007l) . 

Whether such a unified dSph evolutionary scenario might 
conflict with standard CDM models can be tested by compar- 
ing the predicted versus actual number of systems in each of 
the above evolutionary states. As shown by S. R. Majewski 
et al. (in preparation), tidal disruption and A-CDM are only 
forced into confrontation should it be found that too many 
MW dSphs can be well represented by MFL, tidally disrupt- 
ing systems. 
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APPENDIX 

SOME COMMENTS ON THE USE OF KING PROFILES 

Although not particularly pertinent to the results of this paper, the topic of King profiles has often been linked to studies of the 
luminous structure of dSph galaxies. In the spirit of the overall review of dSph tidal disruption models we have presented, a few 
comments about King profile s seem appropriate - 
Hi storically, King profiles ([King 1962, 1966) ha ve frequently bee n used to represent the stellar density distributions of dSphs 
(e.g., Hodgdll966tlEs]mdge1ll988allb1 : lrrwin & Ha tzidimitriou 1995). In this regard it is important to distinguish between [King 
( 1962) profiles, which are a simple empirical description of the distribution of light in a spheroidally-shaped system, and lKing 
( 1966) models, which are physicall y-motivated and derived explicitly for the case of densely populated globular clusters placed 
in an external gravitational field. Kingj (I 19661) showed that under these circumstances, his single-component models evolved 
naturally to the empirical profile he derived earlier for globular cl usters, one whose structure can be described in terms of two 
linear physical scales: a core and a truncation radius (identified bv lKingiri962l as the tidal radius). In globular clusters two-body 
stellar relaxation dominates the internal dynamics of the system and leads to the King profile in an external field, but two-body 
relaxation is much too inefficient in low density dSphs. Nevertheless, King profiles have frequently been found to provide an 
overall adequate empirical parameterization of the structural properties of (the most luminous, inner parts of) the light distribution 
in dSphs. In this usage, the outer physical scale is associated with the asymptotal radius at which the fitted two-parameter profile 
approaches zero (or would approach zero, in the absence of the "break population"; we have identified this empirically -defined 
radius as rr lm , the limiting radius, because its connection to a true tidal radius is undefined (and, indeed, one of the controversial 
questions about dSph structure — e.g., see §4.7). 

In this paper we have used King profiles only in this way — as a handy empirical descriptor simply to define a scale for 
comparing the sizes of the bound central cores of our simulated satellites to each other and to the luminous inner part of Ca- 
rina itself. In this case, a King profile is a logical expediency motivated by the fact that our initial particle configurations are 
Plummer models (i.e., non-tidally-truncated), and as they become shaped by the tidal field they come to resemble King profiles. 
But none of the conclusions in this paper depend on the use of a King profile or any assumptions regarding its validity. With 
the growing popularity of multi-component mass models for dSph systems, simple King models — designed for dense, tidally 
truncated, single-component systems — would naturally be considered inadequa te, and, as a resu l t, the King profile descrip- 
tion has been increasingly r esisted and criticized when applied to dSphs (e.g.. iKleyna et alJl2001b iGilmore et alJl2006l 120071 : 
iKoch et ai1l2007atlWull2007h . Obviously, in these circumstances King models are physically meaningless: If Galactic dSphs are 
not single-component bodies, but rather more complex systems where the dark matter is distributed over much larger physical 
scales than the light, then the physical meaning of r/,„, as a tidal truncation radius in the corresponding King profile is lost — and 
consequently, so too is the interpretation that "breaks" in the density distributions indicate the presence of tidal debris. 

Koch et al.'s (2007a) discussion demonstrates various of the issues at play: "It has frequently been suggested in the literature 
that the detection of member stars well beyond the formal tidal limit radius of dSphs, canonically defined via a single-component 
King model representation of the observed surface brightness profiles, is evidence for physical tides. A large number of dSphs 
have been claimed to show hints of extratidal stars and thus tidal disruption (Irwin & Hatzidimitriou 1995; Martmez-Delgado 
et al. 2001; Palma et al. 2003; Munoz et al. 2005, 2006a, 2006b). Evidence for this has also been reported for Leo I (S06). 
What such studies show, in fact, is that the presumed parametric fit to the measured surface brightness (and/or direct star counts) 
is inappropriate at large radii. Given that there is no astrophysical basis to the application of a King model to a galaxy, an 
astrophysical interpretation of the corresponding labels as physical processes is fraught with peril." 

Unfortunately, the above description oversimplifies the arg uments made to argue in fav or of the tidal disruption hypothesis by 
at least some of the papers they quote. For example, in the iMunoz et alJ (120051 12006a b) papers cited above (and elsewhere), 
additional information — such as the dynamics of the stars, their two-dimensional sky distribution, the implied dark matter 
masses needed to keep the stars bound, and the fact that the radial density profiles resemble models of disrupting dSph galaxies 
— were also weighed into those analyses and provide more powerful and more important support for the notion that these systems 
may be tidally disrupting. In the end, while the observed "breaks" in the density distributions of dSphs are not an unequivocal 
indicator of tidal stripping, they are at least consistent with the presence of tidal debris in the outer parts o f these satellites as 
observed in numerous jV-body simulations, even when multi-compone nt models are used for the dSphs (e.g jRead et al.| [2006b; 
iKlimentowski et al.ll2007l: iPenarrubia. Navarro. & McConnachidl2007l) . 

Finally, it is worth pointing out that there might be circumstances where King profiles (or ot her two-para meter density law) 
are actually physically motivated and therefore are an appropriate physical descriptor of dSphs. I King! d 19621) himself addressed 
this issue by pointing out that his empirically derived density law is applicable to relaxed systems regardless of the relaxation 
mechanism. For example, violent relaxation might provide an alternative to two-body relaxation, perhaps as a process during the 
early formative stages of dSphs. Alternatively, as we have explored in this paper, dSphs might indeed be tidally truncated objects 
where their flat/rising velocity dispersion profiles are the result of tidal stripping, and thus with a truncation radius evolving in the 
luminous matter naturally. Such a situation may be entirely consistent with A-CDM expectations that these systems are (initially) 
embedded in larger dark matter halos, if the systems evolve to where they are reasonably described by single component models 
and King profiles only in the late stages of tidal evolution, a f ter the dark outer layers of these subhalos have already been stripped 
away 21 (see §5.2). Along these lines, IKlimentowski et ail (120071) have shown that a two-component satellite with an initially 
more extended dark matter distribution can naturally evolve through tidal stripping into a system with a radial constant M/L. 

21 It might be argued that, even after significant tidal stripping, the dark matter cocoons surrounding the dSphs' stellar components should keep the initial, 
expected Navarro, Frenk, & W hite (NFW, 1997) density distribution. However, indirec t observations show that the dark matter distribution in dSphs is more 
likely to be cored than cuspy IKleyna et a l. 2003; Goerdt et al. 2006; Gilmore e t alJ2007l) . 



